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The development of quantum walks in the context of quantum computation, as 
^ ' generalisations of random walk techniques, led rapidly to several new quantum 

\^ , algorithms. These all follow unitary quantum evolution, apart from the final 

measurement. Since logical qubits in a quantum computer must be protected from 
decoherence by error correction, there is no need to consider decoherence at the level of 
algorithms. Nonetheless, enlarging the range of quantum dynamics to include 
\^ , non-unitary evolution provides a wider range of possibilities for tuning the properties of 

' quantum walks. For example, small amounts of decoherence in a quantum walk on the 

line can produce more uniform spreading (a top- hat distribution), without losing the 
quantum speed up. This paper reviews the work on decoherence, and more generally on 
non-unitary evolution, in quantum walks and suggests what future questions might prove 
P5 ' interesting to pursue in this area. 
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1. Introduction 



The study of quantum versions of random walks is easily motivated. They have provided 
the engine for several new quantum algorithms, they are of mathematical interest in their 
own right, and in physical systems they form a simple example of coherent quantum con- 
trol over single atoms or photons. They have also found applications from the early days 
of qu antum computation as the clock mechanism in a Feynman computer (jFeynmanl . 
19861) . The role of decoherence requires more explanation. In physical systems we must 
of course consider the effects of decoherence when designing our experiments, but this 
alone does not qualify it as a subject for more general study. Decoherence plays a funda- 
mental role in transforming from the quantum to the classical regime. Quantum walks 
are simple systems in which the intricacies of decoherence can be studied both analyti- 
cally and numerically, making connections with fields such as quantum chaos. Taking a 
broad view of decoherence as any process that tends to reduce quantum coherence, we 
find that it provides a method for tuning quantum random walks to improve their algo- 
rithmic properties. Instead of considering only pure quantum dynamics, we can include 
non-unitary operations in our quantum random walk algorithms, thereby enlarging the 
toolbox for controlling their behaviour. 

The origins of quantum walks can be traced back to t he dynamics of quantu m diffusion, 
well-studied in the physics literature, see for example. iFevnman et al.l (jl964l ). mostly as 
models for physical particles moving on regular lattices. Studies of quantum dynamics 
using a dis crete tirne step described as "quantum r andom wal ks" appeared from the lat e 
1980s fromlCudde^ |l988l ). lGrossing and Zeilingeil |l988h andlAharonoy. Y et al1|l992h . 
the latter with quantum optical applications in mind. Meyeil ('l996b[l st udied similar 
systems with the aim of defining quantum cellular automata (jMeveii . Il996al ) . The current 
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surge of interest in the context of quantum information started with lFarhi and Gutmann 



1998 ) studying continuous-ti me q uantum walks, and continued w ith lAharonov. D et al 



20011 ) . lAmbainis et alj (|200ir ) and lNavak and VishwanathI (j2000l ) studying discrete-time 



quantum walks, all with the goal of applying them to quantum algorithms. As with 
classical random walks, there are a wide rang e of possible dynamics that fit within the 
general concept, for example, Gottliehl (2004); Gottlieb et al. ( 20051 ) defined a discrete- 
time walk with a continuous space and a coin with a continuous set of states. 

In this review we will start with a short overview of purely quantum versions of random 
walks as studied by the quantum information community in their quest for new quantum 
algorithms, then in later sections consider the added effects of decoherence. We will cover 
first discrete-time, coined quantum walks in fj2l then continuous-time quantum walks 
in ^ We mention algorithmic applications in ^ briefly describing two of the known 
algorithms using quantum walks. In ^and ^we consider the effects of decoherence in 
quantum walks, focusing mainly on its use as a tool for enhancing computational speed 
up. There are several proposals for implementing a quantum walk directly in a physical 
system, and most of these studies also contribute analysis of the decoherence relevant to 
the physical system in question: we provide a short summary of these studies in fjT] In 
general, we deal only with quantum walks on undirected graphs, with a brief summary 
of what is known about the largely unexplored realm of directed graphs in 

This review will assume a basic familiarity with quantum mechanics and quantum 
information theory. Those wishing to get up to speed on these areas in the context of 
quantum walks are referred to the excellent and comprehensive introductory review from 



Kempel ( 2003a ). which assumes only a basic knowledge of quantum mechanics. Those 



desiring a gentle overview of quantum walks in the context of quantum computing are re- 
ferred to lKendonI (j2006bf ). A brief but accessible review of quantum walk algori thms from 
a com puter science perspective (no mention of decoherence) may be found in lAmbainis 
(|2003l) . 

The mathematical results and proofs for quantum walks without decoherence have been 
presented in great detail and elegance in the original literature: this review summarises 
what is known without presenting formal proofs but with ample references for those 
interested to pursue them further. For decoherence effects fewer analytic results have 
been derived: those that exist are presented in more detail, along with summaries of 
many related ideas, highlighting the myriad of open avenues that remain to be explored. 



2. Coined (discrete-time) quantum walks 

In this section we describe quantum walks taking place in a discrete space of positions, 
with an evolution using discrete time steps. Just as there are many ways to express the 
dynamics of classical random walks, and many variations on the basic "drunkard's walk" , 
the plethora of different studies of quantum walks have proliferated an equally varied set 
of notations. For this review, we have chosen one of the more commonly used approaches 
and we try to maintain a consistent notation, with mention of alternative methods where 
appropriate. To keep the notation simple, we have taken the liberty of using the same 
symbols for common quantities throughout, with the precise definition implied by the 
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context. For example, P{x, t) is a probability distribution over a set of positions denoted 
by a; at a time t during the evolution. In the section on the quantum walk on the line, 
P{x,t) is for the walk on the line, in the section for the walk on a cycle P{x,t) is for 
the walk on a cycle, and so on. Extra dependent variable may appear where needed: 
P(x, a, t) is a probability distribution with the coin state a also specified. 

We begin with a brief description of one of the simplest examples, the quantum walk 
on the line. This will serve both to set up our notation and to indicate the methods 
of solution. It will also illustrate the first key difference between quantum and classical 
random walks, the enhanced rate of spreading. This is followed by a discussion of quantum 
walks on cycles, illustrating several more key differences between classical and quantum 
random walks. We will finish this section with a brief review of coined quantum walks in 
higher dimensions, on general graphs and regular lattices. 



2.1. Coined quantum walk on an infinite line 

In much the same way as we now know almost everything about the properties and 
possible states of two qubits - though quantum computers will clearly need far more 



ouehlv studied (e.s.. see Ambainis et al.l. 2001: B 


ach et al. , 2004 : Yamasaki et al. 


. 2002: 


Kendon and Treeenna, 


20031 Brun et al.. 


2003f 


a; 


Konno et al.. 20041: iKonno. 2002 


: Carteret et al. 


2003f) . though there is no sueeestion that it will lead to useful quantum walk algorithms 



by itself. 

First, recall the classical random walk dynamics we are looking for a quantum coun- 
terpart of. The walker starts at the origin of an infinite line of discrete points (labelled 
. . . —2, —1, 0, 1, 2 ... ) and tosses an unbiased, two-sided coin. If the coin lands "heads" 
the walker makes one step in the positive direction; if it lands "tails" the walker steps in 
the negative direction. This process is repeated T times, and the position of the walker 
is noted: — T < x < T. If the random walk is repeated many times, the probability dis- 
tribution obtained P{x, T) is binomial, as is well known and easily shown. The standard 
deviation of P{x, T) is VT, i.e., the walker is found on average \/T steps from the origin 
after T steps of the random walk. 

One obvious approach to creating a quantum counterpart of a classical random walk is 
to have the quantum walker follow all possible classical random walks i n superposition . 
This is not feasible in the discrete-time quantum walk, as was shown by Mever ( 1996bh . 



because it is not reversible, and therefore not unitary (which all pure quantum dynamics 
is required to be). We thus make the quantum walk dynamics as similar as possible to the 
classical random walk in a different way, by using a quantum coin. Historically, the role of 
the quan t um co in has caused much confusion, and an equivalent formulation presented by 
Watrouj ( 2001 ) using directed graphs seems to be preferred by some computer scientists. 



For brevity and clarity, in this review we will stick to the formulation using a quantum 
coin, with brief mention of the Watrous variant in SjS] The quantum walker on the line 
starts at the origin and tosses a two-state quantum system (a qubit). The qubit coin - 
also called "chirality" with the two states labelled "left" and "right" in many papers - 
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can be in a superposition of both states, say 

m = 0))=a\0,-l)+(3e''^\0,+l), (1) 

where the "ket" |-) is a standard quantum notation to indicate the complex vector of a 
pure quantum state. On the r.h.s. the kets are basis states: the first entry is the position 
on the hne (a; = in this case) , and the second is the state of the quantum coin, where 
we choose to label the two states |±1). There is an arbitrary phase < 4> < 27r, while 
a and /3 are real numbers, with normalisation + = 1. The quantum walker then 
steps in both directions according to the state of the quantum coin, giving 

|V(t=l)) = a|-l,-l)+/3e^1l,+l). (2) 

Keeping to a pure quantum (i.e., unitary) evolution for the moment, the coin toss and 
step can be written as unitary operators C and S. A single step of the quantum random 
walk is thus 

m+i)} = scm)), (3) 

where can be expanded as a superposition of basis states, 

m)) = Y.''-At)\x,c), (4) 

with complex coefficients ax,c{i) satisfying ^ |aa;.c(t)P — 1. After T steps we have 

|^(r)) = iscfim)- (5) 

The shift operator S has already been specified implicitly by equations ([T]) and ([2]), it 
can be defined by its operation on a basis state \x, c), 

Sjx, c) = \x + c, c). (6) 

It remains to specify the coin toss C, which can in principle be any unitary operator on 
the space of a qubit. The equivalence up to a bias (probability ?y, < 77 < 1, to move in 
the positive direction, (1 — 77) in the negative direction) of all coin operators in the walk 
on a line has been no ted by several authors ( Ambainis et al.l . 2001 : Bach et al. . 2004 : 



Yamasaki et al. . 2OO2I ). Expressing the unitary operator as a matrix we have 

is thus the only possible type of coin for the quantum walk on a line. For rj — 1/2 
(unbiased), is commonly known as a Hadamard operator. 



.(H) _±f ^ 1 

~ V2 U -1 

thus the simple quantum walk on a line is also known as a Hadamard walk. Trivial 
cases 77 = 0,1 give oscillatory motion and uniform motion respectively. The full range 
of behaviour is obtained by choosing different initial coin states, i.e., varying a, /3, (fi in 
equation ([1]). 
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Fig. 1. Comparison of a classical (dots) random walk and a quantum (crosses) 
walk with initial coin state | — 1) and a Hadamard coin operator on a line after 100 
steps. Only even points are shown since odd points are unoccupied at even time 
steps (and vice versa). 



The quantum walk on the line can be solved analytically in many ways. Straightfor- 
ward generalisations of methods that work well for class ical random walks, such as path 
counting, and Fourier transformation were both used by Ambainis et al.l |200l[ ). We will 
give exam ples using thes e meth ods in ^ Path counting (path integrals) was further 



refined in Carteret et al. ( 2003h. and a third me thod using the algebra of the matrix 



operators was presented by Konnol tooi l2005ah Soh ition using the tools of classical 
(wave) optics can be found in Knight et al. ( 20031 2004). Romanelli et al. ( 2003 ) analyse 
the walk on the line by separating the dynamics into Markovian and interference terms, 
allowing an alternative method of solution for the long-time limits, which they relate to 
the dynamics of a kicked rotor. They also show that in the continuum limit one obtains 
the diffusion equation with added interference terms. 

The best way to appreciate the most interesting properties of the solution is in a graph, 
shown in figure [H of the probability distribution of the position of the walker, 



P(x,T) 



+1 

E 

c=-l 



\{cmm 



(9) 



obtained in the standard way by tracing out the coin and taking the square modulus 
of the wavefunction. The quantum walk looks nothing like the classical random walk, 
it spread s out much fas t er in a spiky distribution that is a discrete form of an Airey 
function (jCarteret et al.l . 120031 ) . It is also asymmetric, with the asymmetry determined 
by the initial coin state. A symmetric distribution can be obtained by choos ing the initial 
coin state as either (\ — l ) ±i\+ l))/\/2, or cos(7r/8)| — 1) — sin(7r/8)|-|-l), seelKonno et al 



(2004); iTregenna et al.l |2003h . The moments have been calculated by lAmbainis et al 



200ll ). for asymptotically large times T for a walk starting at the origin, 



a'{T) = {x') = (1 - l/V2)T 
independent of initial coin state, and 

(x)a = a(l - 1/V2)T, 



(10) 



(11) 
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where a € {1, —1} is the initial coin state. The standard deviation (from the origin) cr(T) 
is thus hnear in T, in contrast to \/T for the classical walk. Moreover, unlike the classical 
random walk, the quantum walk evolution depends on the initial state at all subsequent 
times. 



2.2. Coined quantum walk on a N-cycle: mixing times 

The A^-cycle is the Cayley graph of the cyclic group of size N. It is also a line segment 
with periodic boundary conditions applied, so there is little extra work to do to apply the 
solution of the quantum walk on the line to the quantum walk on the cycle. In general, 
when a quantum walk occurs on the Cayley graph of some group, it simplifie s greatly on 



consideration of the Fourier space of the particle (jAharonov. D et al.l . l200lh . On cycles 
we are interested not in how far the walker strays from its starting point, but in the 
mixing properties of the distribution on the cycle. Classical random walks on a cycle mix 
to a uniform distribution in a time O(iV^). That is, after 0{N'^) time steps, the walk has 
all but forgotten its starting state and is equally likely to be found on any site around 
the cycle. More formally, we can choose any small e and find the mixing time M(e) such 
that 

M(e)=min{T | V t > T : \\P{x,t) ~ P^Wtv < (12) 
where P„ is the limiting distribution over the cycle, and the total variational distance is 
defined as 

\\Pix,T) - EE ^ |F(.T,T) ~ P„|. (13) 

X 

In general, any reasonable distance function will do the job here, but we will stick with 
one commonly used in the literature so we can make quantitative comparisons. For a 
classical random walk on a cycle, M(e) ^ 0{N^ log(l/e)). 

In the quantum walk on the cycle, the first observation is discouraging: quantum 
random walks are deterministic, they oscillate forever and in general do not mix even 
instantaneously. But by defining a time-averaged distribution, 

I T 

P{x,T) = -Y^P{x,t), (14) 

t=i 

quantum walks on cycles do mix. Operationally, this just means choosing randomly a 
value of t between 1 and T, then measuring the position of the quantum walker after 
t steps. However, unlike the classical random walk, the limiting distribution lim{T 
oo\{P{x,T)) for a quantum walk is not, in general, the uniform distribution. This is in 
stark contrast with a classical random walk, which always mixes to the uniform distri- 
bution (on a regular undirected graph). 

Exactly as for Af(e) above we can then define the mixing time for P{x,T), which we 
will denote M(e), 

M(e) = min{T | V t > T : ||P(M) - ^«l|tv < e} ■ (15) 
Aharonov. D et al] ( 2001 ) proved that the coined quantum walk on a cycle with an odd 
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number of nodes does mix to the uniform distribution, and has M(e) bounded above by 
0(e~'^A^log A^), almost quadratically faster (in TV) than a classical random walk. They 
also proved a lower bound on the time-averaged mixing times for quantum walks on 
general graphs of bounded degree, suggesting a quadratic improvement over classical 
random walks is the best that can be achieved. 

Notice that there is a price to be paid for the time-averaging: the scaling of the mixing 
time M(e) depends on 1/e, compared t o log(l/e) for a classical r andom walk (where 
no time-averaging of P{x,T) is needed). lAharonov. D et al.l (|200ll ) show that this can 
be avoided by including an amplification step. The quantum walk is run several times, 
each time starting from the final state of the previous walk. Applied in the optimal 
way, the ir bound on th e quantum mixing time reduces to 0(A^log A^log(l/e)). Recent 
work by Richter (j2006al lbl) improves this result to 0(A^ log(l/e)). Since the intermediate 
measurements render the overall quantum walk dynamics non-unitary, we will postpone 
more detailed discussion to §5.51 

This speed up in the rate of mixing of quantum walks is a second key difference be- 
tween quantum and classical random walks. We can also say a little more about their non- 
classical limiting distributions. For example, Hadamard walks on cycl es with an odd num 



berof nodes converge to the uniform distribution (as was proved by lAharonov. D et al 



20011) . but those with an even number converge to a non- uniform distribution unless an 
extra phase S is added to the Hadamard coin operator 



Co = 



-Vv 



(16) 



as 



Tregenna et al.l (|2003^ and iBednarska et al.l (j2003l ) show. By appropriate choice of 



coin operator, a walk on any size cycle can be made to converge either to a uniform or to 
a non-uniform probability distribution. In classical random walks, the properties of the 
limiting distribution depend solely on the form of the graph. 



2.3. Periodicity in coined quantum walks on cycles 



A further curious property of quantum walks on cycles is that for a few special choices 
of cycle size N , a perfectly periodic walk arises where the quantum walk returns exactly 
to its starting state after a fixed number of steps r2, and then repeats the sequence 
over again, returning at 2f2, etc.. The classical random walk has no such behaviour, and 
returns to its initial state only after irregular numbers of steps. This periodicity is not 
connected with whether the limiting distribution is uniform or not, since here we are 
concerned with exact return to the initial state, rather than the time-averaged quantity 
in equation ([T^ . Some of these periodic walks also mix instantaneously, if we allow a 
walk on an even cycle to be considered mixed on just the odd or just the even-numbered 
sites (the same issue arises when considering classical random walks on even-sized cycles) . 

Using a Hadamard coin, the "cycle" of size = 2 is trivially periodic, returning to 
its original state after two steps. A cycle of size N = 4 has a period of ei ght steps. This 
was first noted by Travaglione and MilburnI ( 2002 ). Tregenna et al. ( 2003 ) continued this 
investigation, and found the cycle with N = 8 has a period of 24 steps, but iV = 16 is 
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chaotic and does not return to its initial state exactly even after many thousands of steps. 
These periodic cycles also exhibit instantaneous mixing half way through on their way 
to returning to the initial state, but instantaneous mixing has not been systematically 
studied in discrete walks on the cycle. If the coin is allowed to be biased, a few more 
periodic examples can be found, N = 6 with period 12, and iV = 10 with period 60. With 
judicious choice of phases in the coin operator, = 3 has a period of 12, and N ~ 5 
has a period of 60, clearly related to = 6 and = 10 respectively, but these were the 
only periodic odd-iV cycles found. It is not known wheth er all periodic quantum walks 
on a cycle have been identified, but iTregenna et all |2003h conjecture that there are only 
a finite number of such solutions and that this is nearly all if not all of them. Neither 
periodicity nor instantaneous mixing has been explored on more general graphs beyond 



the tr ivial extensions of the above to the torus and twisted torus in ITregenna et al 



(|2003h 



2.4. Coined quantum walks on higher dimensional graphs 



Classical random walks are not limited to one-dimensional structures, and neither are 
quantum walks. All that is required for the discrete-time quantum walk is a coin that is 
large enough to handle the number of choices at each vertex the quantum walker might 
land on. 

Consider a general graph G, with N vertices in a set V, connected by edges from the 
set E. How one proceeds depends on what prior knowledge one has about V and E. Since 
we are discussing only undirected graphs, if vertices x,y & V are connected by an edge 
Sxy G E, then also ey^ G E and one is allowed to travel both from x to y and from y to 
X. The usual way to represent the structure of the graph is in an adjacency matrix A, 
which has unit entries for each A^y for which Cxy G E, and zeros everywhere else. Since 
Axy — Ayx for an undirected graph, A is symmetric. This representation of the graph 
assumes that Cxy is unique, i.e., there is at most one edge between any two vertices in 
G. Classically, one can subsume multiple edges between the same vertices into a set of 
"edge weights" . However, a quantum walker might traverse both edges at once but with 
different phases, which would in general have a different outcome to one weighted edge. 

Given no further information about G, there could be as many as N{N — l)/2 edges, 
corresponding to the complete graph, for which every vertex has A'^ — 1 edges leading from 
it. In this case the coined quantum walk needs to use a coin of at lea st this si z e (A^ — 1), 
the details of how to implement this quantum walk may be found in iKendonI (|2006al ) . If 
one knows the maximum degree of the graph, d, then one only needs a coin of size d. The 
details of how to imple ment this quantum walk , originally from Watrous ( 20011 ). are more 
accessibly described in Ambainia (2004) and Kendon and SandersI ( 2004 ). A different 
approach, usin g self- loops to make the degree of the graph constant, is mentioned in 
Kempel |2003al) . 



Here we summarise the approach given in iKendon and Sanderd (|2004l ) . First we define 
our Hilbert space for the quantum walk, Tivc which contains the A'^-dimensional Hilbert 
space 

TYv = span{|a;)v : x G ZAr,v(a;|x')v = 5xx'} C T^vc (17) 
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of vertex states, and for the coin a d-dimensional Hilbert space 

Tic = span{|c)c : c £ Z<j} and c(c|c'>c = Sec'} (18) 

where d is the degree of the graph. The basis states of H^c are given by 

Svc = {\x, c) = |a;)v|c)c; a; g Zat, c e Z^} (19) 

with cardinahty Nd. For a basis state |a;, c), the index a; identifies the vertex number and 
c the c*** state of the coin. For an edge e^x' we associate the coin state c with the edge 
at Vx, and the coin state c' with the other end of the edge at Vx'- The values of c and 
c' are arbitrary but fixed throughout the quantum walk, to ensure the quantum walker 
traverses the the graph in a consistent manner. We define the mapping 

C : Zat X Zd ^ Zat X Zd : [x, c) ^ C{x, c) = [x' , c'), (20) 

where {x, c) and {x' , c') label each end of Cxx' ■ 

The unitary quantum walk evolves by repetition of two steps: a coin toss and a condi- 
tional swap. The coin operator 

C : Hvc ^ Hvc : c)^Yl <^"ck: c)c (21) 

is a block diagonal matrix with each block labelled by x. The x-dependence of the coin 
matrix allows sufficient freedom in the quantum walk dynamics for the quantum coher- 
ence properties of the coin to vary between vertices, for vertices to act as origins and 
endpoints, and for vertices to have different degrees from each other. If Vx has degree 
dx < d, we require C^g — for all c values not used to label an edge at Vx ■ This restricts 
the coin operator so it only produces states that have a valid mapping under C. If 

C:, = C^^x,x', (22) 

we have the special case of a fixed degree graph where the coin operator is identical for 
all vertices, as in the walk on a line or cycle described in the previous sections. 
The unitary conditional swap operator is given by 

S:n.c^'H^c-\x,c) ^\x',c') , (23) 

which updates the position of the walker and the coin state according to the mapping 
^ in equation (j20p . i. e., moves the walker and coin to the vertex Vx' along edge Cxx'- 
We note that, by our stipulation that c and c' label opposite ends of Cxx', it follows that 
S = S~^, and is thus unitary as required for quantum evolution. The sequence of a coin 
flip and a conditional swap is a transition over the unit time step, which we denote by 
unitary U = SC. The quantum walk can then be written 

|V^(t))=U*|V;(0)), (24) 

where li^it)) and |'0(O)) can be expressed in basis states as per equation (j4]). 
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2.5. Coined quantum walks on regular lattices 



Analytical solutions of equation p4)l can generall y only be obtained for s pecial graphs of 
fixed degree and hi gh symmetry (see, for examp le, iGrimmett et al.l . 120041 ). The scattering 
theory methods of iFeldman and Hillervl (200J) are worth noting as a possible exception, 
though few applications have been presented. Numerical studies are less constrained, due 
to the simplicity of the step by step evolution of quantum walks. Examples of quantum 
walks on variou s grap h s of degree larger tha, n two can be found in lMackav et al.l ([2002); 



Tregenna et al. 



of Shenvi et al 



(200i); ICarneiro et al.l (|2005i ) and in the quantum walk search algorithm 



|2003 ). see §4.11 Many of the basic results for quantum walks on the 



line and cycle, such as faster spreading, carry over to higher dimensiona l graphs. This 



was s hown numerically for regular lattices of degree three and four by iMackav et al 



ilowever, higher dimensional coi ns have a much wider s et of possible types of coin 



operators, and further numerical study (jTregenna et al.l . 120031 ) shows that the quantum 
speed up is not automatic for all possible choices. 

The role of the lattice symmetry combined with the symmetry in the dynamics of the 
quantum walk i n determ ining the spreading has been further studied and clarified by 
Krovi and Brun |2006al |bl). Most of the examples studied to date have a high degree of 
symmetry, and it makes sense to choose a quantum coin operator that reflects the sym- 
metry in the problem. Two examples are worth noting. First, the DFT (discrete Fourier 
transform) coin operator is unbiased, but asymmetric in that you cannot i nterchange the 
labels on the directions without changing the coin operator. First used by iMackav et al 
(|2002f) . for d = 3 it looks like 




where 



and 



(25) 



are the complex cube roots of unity, For d = 2, the DFT coin 



reduces to the Hadamard coin, equation ([8]), though this is not t he only way to gener- 



alise the Had amard coin to higher dimensions (see, for example, iTregenna et al.l . 12003 : 
Tadei and Zvczkowskiibooel ). Quantum walks using the DFT coin operator have interest- 
ing no n -classical properties th a t have been studi e d bvlMackav et al.l (|2002l ) ; ITregenna et al. 
(l2003l) : ICarneiro et all (|2005l) : iKrovi and BrunI (|2006af ). but none have yet provided any 
quantum algorithms. 

The second com monly used co in operator is a highly symmetric coin based on Grover's 



diffusion operator ( Grover . 1996 ^ with elements 2/d—Sij, shown in matrix form for d = 3, 



AG) 




(26) 



The Grover coin is biased but symmetric, it is the symmet ric unitary oper ator furthest 
from the identity. It was first used in quantum walks b y IWatrousI (120011). a nd is the 
key ingred i ent in the quantum walk searching algorithm of Shenvi et al. ( 20031 ) , see §4.11 
Inui et al.l (|200J) studied the localisation properties related to searching on a two di- 
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mensional lattice. Szegedv ( 2004al fbh introduces a generalisation of the Grover coin that 
quantises an arbitrary Markov chain: essentially this allows for edge weights on the graph, 
and works for graphs of variable degree as well as regular graphs. 



2.6. Coined quantum walk on the hypercube: hitting times 



In one of the few analyt ical studies of coined quantum walks in higher dimensions, 
Moore and Russel] ( 2002h use a Grover coin for the hypercube, which has 2" vertices 
each of dimension d = n. The n-dimensional hypercube is the Cayley graph of , so the 
solution follows the same general me thod as for the line and th e A''-cycle, using Fourier 
transformation of the position space. iMoore and Russell (j2002l ) determine from this so- 
lution that the discrete-time quantum walk on the hypercube of size n has approximate 
instantaneous mixing times. 



Mi„,t(e) = {i : \\P{x,t) - Pu\\tv < e} 



(27) 



for t = nkn/4: for all odd fc > and e = 0(n-^/S). This is an improvement over a classical 
random walk which mixes in time 0(n log n), but requires measurement of the quantum 
walk at exactly the right time (else the walk "u nmixes" again as it procee ds). They also 



consider P{x,t), and find using the methods of lAharonov. D et al.l (|200l[ ). that M(e) is 
exponentially large, 0(2"). 

These results for mixing times are disco uraging , but there is another useful property 
one can test on a hypercube: hitting times. iKempd (|2003bl . 12005! ) proved that a quantum 
walk can travel from one corner of a hypercube to the opposite corner exponentially faster 
than a classical random walk. Kempe defines two different hitting times: instantaneous, 
where one measures the destination corner at the optimal time, and concurrent, where 
one checks the destination after each step to see if the walker arrived yet. This gives us our 
first taste of non-unitary evolution in a quantum walk, since measuring the destination 
no de at each step destroys some of the coherences in the quantum state. Recent work 
bv lKrovi and BrunI (|2006ar ) expands on these ideas, we will describe them in more detail 
in ^5.61 There are other classical algorithms that can cross a hypercube efficiently, so it 
does not provide a quantum algorithm with a true advantage over classical. 



3. Continuous-time quantum walks 



Co ntinuous-time quantum walks on a discrete lattice have their origins back as far 



— .Fevnman et al. ( 19641). Their use for quantum algorithms was first suggested by 
Farhi and Gutmann "^998) , who showed numerically they can reach the ends of certain 



network configurations more efficiently than classical random walks. A proven exponen- 
tial speed up in a quantum a l gorith m using a continuous-time quantum walk came a few 
years later from lchilds et al" (2003), which we will briefly describe in M.2\ In this section 
we will first describe the continuous-time quantum walk dynamics on a general graph, 
then compare it (numerically) with the discrete-time walk on the line. We then give two 
further examples, on cycles and hypercubes, that we will need later. 
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3.1. Continuous-time quantum walks on general graphs 



The co ntinuous-time quantum walk naturally works on any undirected graph. iFarhi and Gutmann 



( 19981 ) simply use the adjacency matrix A, which is symmetric for an undirected graph 



and therefore Hermitian, to form the Hamiltonian for the evolution of the quantum state, 

*l(^I^W)=E(^|H|y>WW>- (28) 

y 

Here H = 7A, where 7 is the hopping rate per edge per unit time, and where is now 

a vector in the Hilbert space of position (vertices in the graph) only (no coin). Continuous- 
time quantum walks achieve what was not possible without a coin in a discrete-time 
quantum walk: the traversing of all possible paths in superposition. The formal solution 
of equation (^5)) is 

|^(i))=e-'^^*|V(0)). (29) 
On the same graph, a continuous-time classical random walk evolves as 

^P(x, t) = 7 E i^-yPiy^ - Ay^cPix, t)} . (30) 
y 

Comparing equation (j28p with equation (j30[) shows that in the quantum walk, the second 
term, necessary for conservation of probability, is missing. We only need a Hermitian 
operator for quantum evolution, and since A^y — Ayx guarantees this, we are free to 
examine this dynamics as well as that obtained in more direct analogy to equation (|30p 
by using the Laplacian, 

d 



i^Mm) = iY.{^^\My){ym)) ~ {y\M^){Am)} ■ (31) 



For graphs where all the vertices are of the same degree d, the Hamiltonian in equa- 
tion ([3T|) becomes H — 7(A — dl), and the solution to this can be written 

=e-*^(^-'^l)*|V(0)). (32) 

Since A commutes with the identity, the two terms in the exponential can be factored, 
giving 

\^{t)) ^e-'i^'e'^'^'^'\il^{0)). (33) 
The factor e'^''^* is only an irrelev ant global phase that makes no difference to observ- 



able quantities if omitted (jAhmadi et al. . 2003i) . thus equation ([33]) is equivalent to equa- 



tion for graphs of fixed degree. For graphs of general degree, however, the dynamics 



with th e second term included will have a different evolution, as IChilds and Goldstone 



(|2004a^ point out. So far, only graphs of fixed degree have been studied in any detail in 
the literature, so the differences between the two versions have not been explored. 

Equation pg]) looks remarkably similar to the discrete case, equation ((24)) : both can 
be written in the form 

|V(t))-U*|V'(0)), (34) 
with unitary operator U = SC for the discrete-time walk and U = e^*^^ for the 
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continuous-time walk. But, unlike the classical case, where the limit of the discrete-time 
walk as the time step goes to zero can be taken in a way that gives the continuous-time 
walk, in the quantum case the similarity is deceptive. The discrete and continuous walks 
have Hilbert spaces of different sizes, since the continuous-time quantum walk has no 
coin: the continuous-time walk therefore cannot be the limit of the di s crete-t ime walk 
unless it can be restricted to a subspace that excludes the coin. StrauchI |2006bl (al) shows 
how to do this for the walk on the line: it turns out the appropriate limit of the dis- 
crete time quantum walk produces two copies of the continuous-time quantum walk, 
one for each of the two coin degrees of freedom. Alternatively, if one adds a coin space 
to the continuous-time quantum walk, a quantum dynamics is obtained that is more 
obvious l y the continuum limit of the d iscrete-time coined quantum walk, see ISzegedv 
( 2004bh : IChilds and Cxoldstonel (|2004bh . The quantised Markov chain formalism intro- 
duced by lSzegedyl (j2004ar ) is yet another convenient wa .y of constru cting quantum walks 
in a parallel manner to classical random walks (e.g., see Weissl . 1994 ). The stochastic ma- 
trix governing the Markov chain corresponds to the adjacency matrix in the graph-based 
description above. 



3.2. Continuous-time quantum walk on the line 

A comparison of the simple case of the walk on the line will show us the main similarities 
and differences between discrete-time and continuous-time quantum walks. Equation (|28p 
has a straightforward solution in terms of Bessel functions, 

+ 00 

\m) = E i-^rum, (35) 

a: — — 00 

where Jx{t) is the Bessel function of order x. Using a hopping rate of 7 = 0.5 so the 
total hopping probability per unit time is one, after 40 units of time we obtain the 
probability distribution shown in figure [21 The discrete-time quantum walk evolved for 
55 time steps is also shown for comparison. We can see that the shapes of the distribution 
are comparable, the difference in the height of the peaks is due largely to the continuous- 
time walk having support on both odd and even sites, while the discrete-time walk is 
restricted to sites with the same parity as the time step. Nonetheless, the continuous- 
time walk has a very oscillatory nature about it, especially in the central region. The 
two types of walks propagate at somewhat different speeds, as evidenced by the different 
time instants at which they have the same width. Both spread linearly, just with different 
constant prefactors, so both show a quadratic speed up in their spreading compared to 
classical random walks. For the continuous-time quantum walk starting at the origin, for 
asymptotically large times T, 

a\T) = {x')=T'^/j. (36) 

The continuous-time quantum walk on the line is always symmetric, if the hopping rate 
is the same for both directions, this is the main difference with the discrete-time walk, 
where the coin generally skews the walk unless special choices are made for the initial 
state. 
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-• cont-time, T = 40 
X discr-time, T = 55 



20 20 

position, X 




Fig. 2. Comparison of continuous-time (dots) and discrete-time (crosses) quantum 
walks on a line starting at the origin after 40 and 55 steps respectively. The 
discrete-time walk has initial coin state (| — 1) -|- i\l))/^/2, and only even points 
are shown since odd points are unoccupied at even time steps (and vice versa) . All 
points are shown for the continuous-time walk, which has hopping rate 7 = 0.5. 



3.3. Continuous-time quantum walk on a N-cycle and other circulant graphs 

The continuous time quantum walk on a TV-cycle is straightforward to analyse. The adja- 
cency matrix A of the iV- vertex cycle graph is a circulant matrix, it has eigenvalues A^; = 
2 cos(27ra;/A'') with corresponding eigenvectors \bx), where {y\hx) = exp{—2Trixy / N) , 
for a; = 0, 1, . . . , iV — 1. Taking the initial state of the quantum walk to be \ip(0)) — |0), 
then \ip{t)) — e^**^|0) can be solved by decomposing |0) in terms of the eigenvectors 
\bx), giving 



x=0 

One can apply the same notions of mixing times for a continuous-time quantum walk as 
for the discrete-time walk. Like the discrete-time quantum walk on a cycle, the probability 
distribution P{x,t) of the continuous-time walk does not mix asymptotically, and is 
known to have exact i nstan taneous mixing only for a few special cases, = 3 and 
A^ = 4 lAhmadi et al.l . l2003h . i.e. even less cases than are known for the discrete-time 



walk on the cycle as mentioned in H2.3I We can define the continuous-time version of 
equation (fT4|) . the average probability 



1 



T 



P{x,T) = - dtP{x,t) (38) 

where P{x,t) — \{x\tlj{t))\'^ . While this does mix asymptotically, like the discrete-time 
walk, in general the limiting distribution is not uniform: it retains one or more peaks 
reflecting the initial state. Notice that the properties of continuous-time walks on cycles 
are not dependent on whether the cycle has an odd or even number of nodes, unlike the 
di screte-time quantum w alk. 

Adamczak et al. I (|2003l) show that continuous-time walks on cycles are nearly uniform 



mixing, i.e., if one relaxes the condition on how well they they approach the uniform 
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distribution from an arbitrarily small e, they satisfy 



M(l/4iV) = min<{T | V OT: \\P{x,t) - 1/N\\ty < —} ^ 0{N). (39) 



Adamczak et al.l (|2003l ) extend their analysis of mixing in continuous-tim e quantum walks 
to a wider variety of circulant graphs, this is continued in lLo et"aL ( 2006 ): r nixing to non- 
unifo rm distributions turns out to be the norm for these types of graphs. ^ Carlson et al 



identify graphs where, with the use of edge weights, any distribution on the graph 
can be obtained (universal mixing). They also identify another class of graphs that are 
instantaneous uniform mixing, the claw or star graphs. 



3.4. Continuous-time quantum walk on the hypercube 

Our next example o f a co ntinuous-time quantum walk is on a n-dimensional hypercube. 



Moore and Russelll (|2002l ) analytically solved the continuous-time quantum walks on the 
hypercube as well as the discrete-time walk (see The analysis for the continuous- 

time walk makes use of the hypercube's product graph structure. We label the vertices 
with n-bit strings, with edges connecting those pairs of vertices that differ in exactly one 
bit. Then, using the Pauli matrix, which is the bit flip operator, 

the adjacency matrix can be decomposed into the sum 

n 

A==^l(g)---(g)Cr:r«)---«)l , (41) 

J = l 

where the jth term in the sum has as the jth factor in the tensor product. Each 
term thus flips the bit in the vertex label corresponding to traversing the edge to the 
appropriately labelled neighbouring vertex. Using H — 7A for the quantum walk, we 
hav^ 

n 

U* = e-'^' = Y[ 1(g) ■■■(g) e-''"'"'' (g) ■■■(g) l=[e-''^"'^]'^" 



cos(7i) —i sin(7i) 
—i sin(7i) cos(7i) 



(42) 



Applying U* to the initial state |?A(0)) = |0)®", we have 

U*|V'(0)) = [cos(7i)|0) - i sin(7t)|l)]®" (43) 

which corresponds to a uniform state exactly when ^t is an odd multiple of j. The 
choice of 7 is thus pivotal for making a fair comparison of the instantaneous mixing time. 

t iMoore and Russelll 1 I2OO2I) use U = e'^, this is just the complex conjugate and makes no difference to 
observable quantities. 
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Moore and Russell choose 7 = so that the total probability of making a hop 



to any neighbouring site per unit time is unity. With this choice, the continuous-time 
quantum walk mixes exactly instantaneously in time 0{n), a logarithmic improvement 
over 0{n log n) for the mixing time of a classical random walk. They also showed that for 
the time-averaged mixing time given by equation ()15|) . the continuous-time walk never 
mixes. These results can be compared with the discrete-time walk on the hypercube, 
which has approximate instantaneous mixing (instead of exact), and M(e) mixing in 
exponential time. 

The continuous-time quantum walk on the hypercube also hits the opposite corner in 
linear time, exponentially faster than a classical random walk: we will obtain this result 
in ^6. 31 as a special case of the calculation of the effects of decoherence. As with the coined 
quantum walk , this e xponentially fast hitting time is highly d ependent on the s ymmetry- 
Keating et al. ( 2006 ) apply results from localisation theory ( Andersonl . 1958 ) to argue 



that this behaviour is exceptional, and the norm on less regular graphs is a quantum 
walk that tends to stay near its starting state. 



4. Algorithms using quantum walks 



Starting with Shor's algorithm for factoring large numbers (jShoii 119971 ). many of the 
quantum algorithms found so far belong to t he same f amily , based on the use of Fourier 
transforms to identify a hidden subgroup (jLomontl . 12004 . provides a recent review). 
This works well (can be exponentially better than known classical methods) for Abelian 
groups, but extending the method to non-Abelian groups, where some of the notorious 
hard problems, such as graph isomorphism, lie, is proving tricky. An obvious place to 
look for new ideas is where classical algorithms are having the most success, to see if 
a quantum version could be even faster. Randomised algorithms are one such arena, 
providing the best known algorithms for approximating the permanent of a matrix 
(IJerrum et al. . 2001), fi nding satisfying assignments to Boolean e xpressions (fcSAT with 
k > 2) (Schoning, 19991). estimating the volume of a convex body ( Dver et al. . 1991 ). and 
graph connectivity ( Motwani and Raghavan . 19951) . Classical random walks also under- 
pin many standard methods in computational physics, such as Monte Carlo simulations. 
Expanding the repertoire of methods for quantum algorithms was the motivation be- 
hind the recent upsurge of inte rest in quantum wal ks. T he first prope r algor ithms using 
quantum walks appeared from Childs et al ( 20o3l and Shenvi et al. ( 2003 ). and more 
have since followed: for a short survey, see Ambainis ( 20031 ). We will briefiv describe the 
first two quantum walk algorithms, since they are of two distinct generic types, and later 
algorithms are mostly variants of the quantum walk search outlined in the next section. 



4.1. Quantum walk searching 

Quadratically faster spreading, as described in §2.1) isn't a quantum algorithm, but 
quantum searching of an unsorted database (e.g., finding the name corresponding to 
a given number by searching a telephone directory) is the reverse process: start in a 
uniform superposition over the whole database and home in on the marked item. The 
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first quantum algorithm for this problem is due to iGroverj (jl996l ) , using a method of 
amplitude amplification, from which a quadratic speed up over classical searching is 
obtained. It c an be shown that a q uadratic speed up is the best possible improvement for 
this problem (jBennett et al.l . ll997l ). A classical search of an unsorted database potentially 
has to check all N entries in the database, and on average has to check at least half. A 
quantum search only need s to make O ( yfN) queries, though the queries ask for many 
answers in superposition. IShenvi et al.l (|2003l ) showed that a quantum walk can also 
search an unsorted database with a quadratic speed up. They represent the database by 
the vertices of a graph with a regular structure (lattice, hypercube...) and start with the 
quantum walker in a superposition of all positions on the vertices. The quantum walk 
proceeds using a Grover coin operator of appropriate dimension at every vertex except 
the marked item. For the marked item, almost any other coin operator can be used, this 
is enough to break the symmetry of the quantum walk and cause the walker to converge 
on the marked vertex in ^ 7r/2-\/iV/2 steps. 

Since then, several variations on quantum searching for a single item have been anal 



ysed, all essentially searching for a set of rnarked items of some sort (jMagniez et al- 
2OO5I : Childs and Eisenbere . 20051 : Ambainil 2004). Spatial searching, where it is also 



counted as a cost to move from one item in the database to the next, is also faster with 
a quantum algorithm, and curiously, the continuous-time walk fi nds this a little harder 
than t he discrete-time walk: in low dimensions it needs a coin too ( Childs and Goldstonei 
2004allbh . 



4.2. "Glued trees" algorithm 



Childs et al.l (j2003[ ) proved that a continuous-time quantum walk could perform expo- 
nentially faster than any classical algorithm when finding a route across a particular sort 
of network, shown in figure [S] This is a rather artificial problem, but proves in principle 
that quantum walks are a powerful tool. The task is to find your way from the entrance 
node to the exit node, treating the rest of the network like a maze where you cannot see 
the other nodes from where you stand, only the choice of paths. It is easy to tell which 
way is forward until you reach the random joins in the centre. After this, any classical 
attempt to pick out the correct path to the exit get lost in the central region and takes 
exponentially long, on average, to find the way out. A quantum walk, on the other hand, 
travels through all the paths in superposition, and the quantum interference between 
different paths allows the quantum walker to figure out which way is forward right up 
to the exit, which it finds in time proportional to the width of the network. The proof 
is quite technical, involving oracles and consistent colourings of the network, and the 
ability to simulate a continuous-time quantum walk o n a discrete quantu m computer. 



We refer the reader to the original paper for the details (jChilds et al.l . l2003h , since we do 
not need them here for our discussions of decoherence. 

Discrete and continuous-time quantum walks are generally expected to have the same 
computational power. They give broadly similar results for algorithmic properties such 
as mixing times and hitting times in those cases where both for ms have been applied to 
the same problem, with the possible exception of spatial search ^Childs and Goldstonei 
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Entrance 




Exit 



12 3 4 



5 6 78 9 



Fig. 3. "Glued trees" graph used in the algorithm of lChilds et alJ (j2003l ). Example 
shown is for = 4, with 2N + 2=10 colmims labelled at the bottom of the figure, 
and 2(2^^+^^ — 1) = 62 nodes. The task is to travel from entrance to exit without 
getting lost in the randomly joined middle section of the graph. The gap between 
columns 4 and 5 is for clarity in the figure and is not significant in the algorithm. 



2004a|), which, as noted above, requires a coin for maximum efficien cy even in the 
continuous-time walk in low dimensions (jChilds and Goldstond . l2004bl ). For a discus- 
sion of how both discrete and continuous-time quantum walk implementations w ould use 
roughly the same computational resources for the "glued trees" problem, see iKendon 
( 2006ah . 



5. Decoherence in coined quantum walks 

As we have seen, while there are many similarities between discrete and continuous- time 
quantum walks, there are also some distinctly different behaviours, and this turns out to 
be true also under the influence of decoherence. We will thus devote a separate section to 
each type, starting here with coined quantum walks and following in ^with continuous- 
time quantum walks. We are going to take a very broad view of decoherence in this 
review as any dynamics that tends to remove the quantum coherences in some way, be it 
unwanted (as in environmental decoherence in an experimental system), intentional (to 
tunc the properties of the quantum walk) , or a byproduct of some other operation, such 
as measur ement. One of the earliest uses of non-unitary quantum walks is measurement- 
based, by Aharonov. Y et al. ( 19921 ). 



Quantum walks are a very br oad class of quantum dynami cs, with overlap into related 



areas such as quantum graphs (jKottos and Smilanskvl . Il997f l. One way to justify a par 



ticular quantum dynamics as being a "quantum walk" is to see if it turns into a classical 
random walk when decohered. Since classical random walks are also a very broad class of 
dynamics, this gets reasonably close to a workable definition of a quantum walk, at least 
for the discrete-time case. Many early studies of quantum walks took the trouble to show 
numerically that, for specific cases, their quantum walks decohered into classical random 
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walks, (e.g., see iMackav et al.l 120021 ) who dephase the coin in their study of a quantum 



walk on the line. A more systematic treatment of the quantum to classical transition 



in a general quantum walk appears in Kendon and Sanders! 1 2004 ). emphasising the im 



portance of demonstrating that quantum walks exhibit both wave (pure quantum) and 
particle (decohered) dynamics and, for a non-unitary quantum walk, being able to inter- 
polate between th ese two modes of behaviour by turning the decoherence up or down. 
Kosik et all ( 20061 ) explicitly calculate the cases of Grover and DFT coins on a Cartesian 



lattice (Z"^), with random phase shifts applied to the coin to induce decoherence, rather 
than measurements. They show that the resulting distributions match the expected clas- 
sical random walks (which may be biased if the initial state or coin operator is biased, 
compare §2.ip . 

Decoherence is usually modelled as a non-unitary evolution of the quantum walk, so we 
will need some extra formalism for mixed states using density matrix operators. Following 
the notation in the (time-dependent) density operator 

P = ^^Pxc,x'c'\x,c){x',c'\, (44) 
x,c x' ,c' 

is a positive (p = with positive real spectrum), unit-trace, bounded linear operator on 
Ti-vc, in the basis Bvc, equation In the decomposition of p into the computational 
basis, the rows and columns of p are indexed by xc and x'c', which run over the position 
states and coin states of the Hilbert space. The state is pure iff p^ ~ p. A typical initial 
condition is p{t = 0) = |0,0)(0, Q| corresponding to the walker starting at vertex vq 
carrying a coin in the state labelled zero. In general the density operator is mapped to a 
new density operator via a completely positive (CP) map 

U : p^Up (45) 

The CP map hi performs both the coin flip and the conditional swap over one time step. 
More explicitly, we can write 

iee iee 

with j an index of non-unitary evolutionary 'instances' and Uj the corresponding Kraus 
operator. These instances may be discerned by a measurement record, with j the record 
index. The cardinality of can be finite, countably infinite, or even uncountable. In the 
case of unitary evolution, 8 has a cardinality of one, so there is a single, unitary U for 
which p 1-^ Up = U pU^ . Unitary quantum walk evolution can be expressed as 

p{t) = U'p{0), U EE 5C, SCp = SCpC^S'' , (47) 

where, for the discrete-time walk, we assume i G Z. Thus, for the unitary walk, a single 
step is given hy U = SC, while for the non-unitary walk we can either add an extra 
non- unitary operation (measurement) {Pj}, 

U-.p^^Fj SCpC^'S^'V] , (48) 

iee 
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or, replace the coin and/or the shift operators by non-unitary operators, depending on 
what sort of decoherence or measurements we are considering. 

Most of this section consists of unpacking equation (|48|) into specific cases, and analysing 
the effects. One may ask, for example, how quickly the quantum walk becomes classi- 
cal as the decoherence is increased. Are quantum walks sensitive to small amounts of 
decoherence, or are the quantum effects robust under environmental disturbance? The 
first studies of decoherence, beyond merely verifying the classical limit is a classical ran- 
do m walk, were analysis of the lik e ly errors in pr o posed experimental imp lementations 

too. 



Travaghone and Milbur J ^200^ : ISanders et aD l|2003t ): loiir et all ^oi). These, 



are focused on the properties of the pure quantum walk, we provide a short summary 
in JJTI The firs t consideration of the algorith mic properties of partly decohered quantum 
walks were by Kendon and Tregenna ( 2003h . who used numerical simulation on a variety 
of discrete coined quantum walks. We will briefly review their results for the walk on the 
line, to provide an overview of the typical effects of decoherence, then examine some of 
the analytical calculations that confirm and expand these initial observations. An alter- 
native approach is to model the entire system of quantum walker plus environment, we 
include examples of this in ^5.4[ where the different methods of decohering the coin are 
compared. 



5.1. Effects in the walk on the line 



First we unpack the superoperator notation of equation (I48p into the specific case of 
decoherence events or measurements occurring independently at each time step. 



p{t+l) = (1 -p)SCp(i)C+S^ +p^PjSCp(t)Ctstl 



(49) 



Here Pj is a projection that represents the action of the decoherence and p is the prob- 
ability of a decoherence event happening per time step, or, completely equivalent math- 
ematically, to a weak coupling between the quantum walk system and some Markovian 
environment with coupling strength p. Equation (|49p lends itself readily to numerical 
simulation since p, S and C can be manipulated as compl ex matrices, while the P,- gen- 
erally remove some or all of the off-diagonal entries in p. iKendon and Tregenna (l2003h 
took equation ([49]) and evolved it numerically for various choices of Pj, projection onto 
the position space, projection into the coin space in the preferred basis | ± 1), and projec- 
tion of both coin and position. Motivated by the likely form of experimental errors, they 
also modelled an imperfect Hadamard by applying a Gaussian spread of st andard devia- 
tion ^/ p7r/4 about the perfe ct value of 7r/2 implicit equation ([8]), compare 

i.r " 



Mackav et al 



(12002) . IShapira et all (|2003h ater modelled imperfect quantum walk op erations such as 



this in more detail, still numerically, and drew the same conclusions, while Konnol ( 2005bl ) 
treated the general case of a randomised coin, showing analytically that the cl assical ran- 



dom w alk is obtained. An imperfect shift by the walk er has been stud i ed in iDiir et al 



|2002h . see gZl and also ( in the form of broken links) bv lRomanelli et al.l (|2004f l. see 

To quantify the change in behaviour in the walk on the line, we can calculate the 
standard deviation, equation (jlOp . which is now also a function of the decoherence rate. 
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Fig. 4. Standard deviation a{T,p) of the p osition on a line for different models of 
decoherence, for T = 100 time steps, (from Kendon and Tregenna . 2003h . 




-20 20 

position (x) 



Fig. 5. A quantum walk on a line of 100 steps is progressively decohered by 
random measurements with probability p per time step as given in the key. for 
p = 0.03, an approximately 'top hat' distribution is obtained. 



In each of these cases iKendon and Tregennal (|2003[ ) found the same general form for the 
decay of a(T,p) from the quantum to the classical value, with small differences in the 
rates, as shown in figure ID The slope of a{T,p) is finite as p — s- and zero at p = 1. The 
decay of the spreading rate (as quantified by a{T,p)), is much as one might have expected 
given quantum states are generally fragile in the face of environmental disturbance. The 
interesting feature is seen in the shape of the distribution of the position as decoherence 
begins to take effect. The changing shape as decoherence is increased is shown in figureEl 
Notice the good approximation to a top-hat distribution between ±T/^/2 that appears 
for p = 0.03. For computational physicists who use random walks to sample distributions, 
this is a very desirable feature, since it provides uniform sampling over a specific range. 
Moreover, this result is only obtained when decohering the position, compare the three 
examples in figured The optimal decoherence ratep„ can be obtained by calculating the 
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Fig. 6. Distribution of the position for a quantum walk on a line after T — 200 
time steps. Pure quantum (dotted), fully classical (dashed), and decoherence at 
rate shown on part of system indicated by key (solid). U niform distribution 
betw een -T/V2 <x < T/V^ ( crosses) also shown, (from lKendon and Tregennal 
20031) . 



total variational distance, equation (|13p . between the actual distribution and an ideal top 
hat distribution. The optimum decoherence rate depends on the number of steps in the 
walk, p„T ~ 2.6 for decoherence on both coi n and position and PuT 5 f or decoherence 
on the position only (obtained numerically in Kendon and Tregenna ( 2003h ). Intriguingly, 
recent work by Malover and Kendon ( 2006[ ) shows that the optimal top- hat occurs for 
the decoherence rate that just removes all the quantum correlations by the end of the 
quantum walk. 

Decoherence on the coin only does not produce a top hat distribution, in stead the distri- 
bution develops a cusp as it passes from quantum to classical. The work of lLopez and 
( 20031 ) provides some insight into w hy decoherence on the coin and on the position pro- 
duce qualitatively different results. iLopez and present an elegant analysis of deco- 
herence in the quantum walk using a discrete Wigner function to bring out the pattern of 
correlations and their decay when the coi n is subject to deco herence. Their decoherence 
model follows much the same approach aslBrun et alj (|2003al ). which will be described in 
§5.21 In their study, Lopez and Paz ( 2003 ^ restrict the dynamics to cycles. However, most 
of their results are for a short numbers of time steps such that the walk has not joined up 
round the cycle, and the results are identical to the walk on the line. Wigner functions 
represent quantum dynamics in a quantum phase space, thus showing the behaviour of 
both position and momentum degrees of freedom in the same picture. This brings out 
beautifully that while decohering the coin state is sufficient to reduce a quantum walk 
on a cycle to a classical walk, if the walker started in a superposition of two position 
states, this remains untouched by the decoherence. Such a walker with a decohered coin 
performs a superposition of two classical random walks, each starting at a different po- 



V. Kendon 



24 



sition. The converse does not apply, because the position is conditioned on the state of 
the coin, equation ([2]), thus decohering the position decohers the coin with it. 

Our first example of decoherence has given us a good guide for what to expect in 
general: besides the rapid degradation of quantum behaviour, there are interesting effects 
in particular ranges of low decoherence. These effects are not algorithmically significant, 
that is, they don't alter the scaling of quantum processing, but nevertheless may be 
useful in practice, and represent intriguing and complex behaviour worthy of study in 
their own right. For practical purposes (assuming that one day we have a quantum 
computer available), the improved top hat profile might still be a useful optimisation to 
get the most out of the computational resources. 



5.2. Dephasing the coin in the walk on the line 

In the next two sections we will consider the analytical treatment of decoherence in the 
walk on the line in some detail. First we consider decoherence only on the coin. This will 
also provide an example of how to use Fouri er transforms to simp lify the dynamics and 



obtain the solution. We follow the method of lBrun et al.l (|2003al |bl) to solve equation (|48l) 
for the coined quantum walk on a line with dephasing applied to the coin. We thus have 
S given by equation ([6]), 

S|a;, c) — \x + c, c), 
where \x) and |c) are basis states, and C by equation ([8]) 

p(H) _ 1 1 M 

Transforming from the position basis |a;) to the Fourier basis |fc) such that, 

E 

we find that 



fc,c)= Ve''=-|x,c) (50) 



Ufe|fc,c) = C,|fc,c) (51) 
where Ck acts only on the coin degrees of freedom. 

We now move to density matrix formulation, for a walk starting at the origin in coin 
state |-0o) (not necessarily a basis state), 

/dk f dk' 

where \k,ipo) is the tensor product \k) (g) lipo)- Using equation pS)) to introduce a deco- 
herence operator, we can thus write 

/dk f dk' 
_ / _|fc)(fc'|^7'L'IV'o)(V'o|, (54) 
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where 

Vkk'\i'o){i'o\ = ^P,Cfe|^o)(V^o|Ctp], (55) 
j 

and the set of projectors {Pj} is the decoherence acting on the coin that we will need to 
specify in order to carry out explicit evaluation of equation ([54|) . 

While it will in general be difficult to evaluate p{t) explicitly, the second moment 
will give us the main information we need to determine how fast the walk spreads. 
Calculating the moments allows considerable simplification even before specifying the 
form of decoherence, 



- j —{x\k){k'\x)Tr, [7'L'l^o)(^o|] , (56) 



where rrc[.] traces over the coin degrees of freedom to aggregate the probability over the 
different coin states. From equation (|50p . {x\k) — e**^^, and we may carry out the sum 
over x, giving 

= JdkJ dk'S^^^ [k - k')Tr, [Vl, l^o) (V'o I] , (57) 

where S^"^\k— k') is the mth derivative of the delta function. This can then be integrated 
b y parts. 

Brun et all |2003bh chose pure dephasing for the form of the decoherence on the coin, 



the coin projectors in equation (|55p are thus 



1 / e±*^ 



(58) 



This can be solved explicitly, giving in the large t limit, 

(x^) - {xf ~ t{cot^ 29 + csc^ 29) + 0(1), (59) 

where 9 measures the strength of the dephasing, 6 = tt/A being complete dephasing 
(classical walk) and = is a pure quantum walk (the approximation above is not valid 
for exactly = 0). This shows that even a small amount of decoherence renders the 
quantum walk classical in the sense that the standard deviation of the position scales as 
y/t, but with a larger prefactor the smaller the dephasing, cot^ 26*^ 1/9^ for small 9. 

We will return to decohering the coin in §5.41 where a variety of other ways to deco- 
her the coin will be considered and compared. Next, we will complete the picture for 
decohering a walk on the line with a calculation of the effects of decoherence in the 
position. 



5.3. Decohering the position in the walk on the line 

When decoherence acts on the position degrees of freedom, the simplification offered 
by the factorisation of the dynamics in Fourier space no longer helps: we would have 
to apply the Fourier transform to the decoherence superoperator too. The following 
calculations of the asymptotic behaviour for small p <C 1 illustrate how real space (path 
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counting) methods may be used instead. Following iKendon and Tregennal |2002l) , we 
calculate a{T,p) analytically for pT ^ 1 and T ^ 1 for the case where the {Pj} are the 
projectors onto the preferred basis {|a,2:)}, i.e., decoherence affecting both walker and 
coin simultaneously. Since, as already explained, decohering the position also decohcrs the 
coin, it makes only a small difference whether or not we also decoher the coin explicitly. 

The probability distribution for finding the walker in the state \a,x) in the presence 
of decoherence can be written, 

P(x, a, r,p) = (1 - pfP{x, a, T) + p{l - pf-'P^'^ (x, a, T) + . . . , (60) 

where P{x,a,T) is the distribution obtained for a perfect walk and P^^^ {x, a,T) is the 
sum of all the ways to have exactly i noise events For example, we can write 

T 

P(i) {x, a,T)^Y.T.T. P^y^ ^' t)Pvb{x, a,T-t), (61) 

t=l y b 

where Pyb{x, a,T ~ t) is the distribution obtained from a perfect walk starting in state 
\y, h) for T — t steps. For the ideal walk, a-^(r) = Yl,x x^P{x, a, T), and for the walk 
with decoherence, 

a2(r,p) = ^^a:2p(x,a,T,p). (62) 

X a 

Taking equation (j60p to first order in p, and substituting along with equation (|6ip into 
equation (p^ gives to first order in pT 

<j^T,p) ~ ^ I (1 - pT)Pix, a,T)+pJ2Yl P^V. b, t)Pyb{x, a, T - t) I . (63) 

x,a I t—1 y^b \ 

The first term on the r.h.s. is (by definition) (1 - pT)a^{T). Noting that Pyb{x, a,T-t) 
is a translation of a walk starting at the origin, Pyb(x, a,T — t) — Pob{x — y,a,T — t), 
and relabelling the summed variable x to {x + y) then enables the sums over x and a to 
be performed in the second term, 



E E p^y^ E E(-^ + vfp^^ix, a,T-t) 

t—l y b X a 

T 

= pEEE^(y'^'^)H(r-i) + 2y(x)|,rny2}. (64) 

t=l y b 

From equation (jlOp. <yQ^,{T — t) does not depend on 6, so the summation over y and h 
may be performed trivially. The remaining summation over y and b applied to y^ gives 
cr^it) by definition. This leaves only the evaluation of 

T T 

^^'E E p^y^ i)y(^)Z''^ = - y^)T.(p -t)Y.y^ p(y^ (^s) 

t = l y,b t=l y,b 

where we have used equation pT|) for {x). We note that this term does not depend 
on whether the initial coin state is plus or minus one and so we may include both 
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these possibilities equally. Also, by the symmetry of the walk, it is possible to rewrite a 
probability function for travelling from state |0, c) to \y, b) in the reverse order, i.e. as a 
probability for moving from \y,b) to |0,c). Care must be taken to ensure that the signs 
of each term due to the coefficient yb in the summation arc maintained. We obtain 

T 

2p^P(y, 6, t)y{x)^^,-'^ = p(l - 1/V2) ^(T - t)J2 (1 - 2Sb,c)yb Pyb{0, c, t) 

t,y,b t—1 y,b,c 



:p{l-l/V2)^(T-i) 



yb PyAo, c,t)-2j2yb Pyb{o, b, t)] 

y,b,c y,b 



(66) 



treating the two parts with and without a delta function independently. Expanding the 
summations over b and translating the position basis by —y gives 

T 

t=l y,b 

T 



= p(l-l/V2)^(T-i) 



E y^o,-i(2/, c, t) - 4 ^ y Pv,+i{y, +1, t) 



(67) 



The final summation over y may be bounded above by noting that 

Po-i(y,-i,<) = |(y,-i|C/*|0,-i)|2 

^ i|(y + 1, l|C/*-i|0, -1) - (y + 1, -l|C/*-i|0, 

c 

Using this in equation ([67|) gives, 

T 

2p^P(y,6,t)y(4^*' < p(l-l/^/2)^(^-^)[(y)U-2(y)(*:^ 



(68) 



t,y,b 



T 



< p{l-l/V2)J2iT-t) (1 - 1/V2)t + 72 



(69) 



where equation pT|) has been used for the average values. Note that equation (fTTj) is 
exact only for asymptotically large times, so the second step in equation (j69p introduces 
further approximations from the contributions to the sum at small times t. Combining 
these results in the full expression for (T, p) , equations ((63l [64]) , and performing the 
summations over t using J2 * = T{T + l)/2 and E = + + T/6, gives 



l-^pr+p(V2-l) 



Taking the square root gives as an upper bound on the standard deviation, 



a{T,p) < a{T) 



pT 



+ ^(l-l/V2) + 0(p2,l/T) 



6\/2 %/2 



(70) 



(71) 
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Equation ((7T|) compares well with simulation data in lKendon and Tregenna ( 2003 ). once 
a second order correction for <t{T) = (1 — 1/V2)^^^(T' — 1/T) is taken into account. 
The bounding procedure applied here is reasonably accurate, numerical studies give the 
coefficient of p in the above expansion as 0.09566, compared with the bound of 0.20711. 
The first order dependence is thus proportional to pT, the number of decoherence events 
during the whole quantum walk. For a given decoherence rate p, the standard deviation 
initial ly decreases linearly in T. This calculation first appeared in lKendon and Tregenna 
(|2002h . 

A similar expansion at the classical end of the full sum in equation (|60p gives 



P{x, a, T, q) = il- qfP^^\x, a, T) + q{l - qf-^P^^-^\x, a, T) + . . . , (72) 

where we have defined q = {1 — p) as the small parameter. We know that a, T) 

is the classical walk, so 

4(T)=r = 5:5:x^p(^)(x,a,T), (73) 

X a 

and we can write 

a2(r,g) = ^^a;2p(a:,a,r,g) 

X a 

= (l-g)V2(T) + g(l-<z)^-i^^x^p(^-i)(x,a,T) + .... (74) 

X a 

However, to see a difference between the classical and quantum walks, we need to have 
four consecutive steps of the quantum walk, since the position distributions are identical 
for the first three steps. The first term that differs from classical is thus derived from 
part of (x, a, T) given by 

P(4=«) (x, a,T) = Y.llll P^^^ ^' ^Kvl) ^' 4)^i+l,c(^, a, T - 4 - t), (75) 

t—1 y,z h,c 

i.e., t classical steps, 4 quantum steps, T — 4 — i classical steps. (Also included in 
, T) are combinations with the 4 quantum steps not adjacent to each other.) 
Since classical segments of the walk arc not influenced by the initial state, we may esti- 
mate the variance of the total contribution of P'^^'^'^\x, a,T) by summing the variances 
of the segments: 

T-4 

'^(4c,)(r) = E [^c W + ^'(4) + <jUT~t-i)]={T-4)[T-4 + a'{i)]. (76) 

t=i 

Explicit calculation (P(a;,a,4) has sixteen terms) gives (T^(4) — 5. From the continued 
expansion of equation (|72p the prefactor is g*(l — q)'^~^, so for small q and large T, we 
have cr^(r, q) ~ r(l + g*) and for the standard deviation 

a{T,q)^VT(l + ^), (77) 



2 



to lowest order in g = 1 — p. 
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A si mpler model of decoherence in both position and coin space is solved bv lRomanelli et al 



(|2004l ). They apply measurements at regular intervals and project the coin into the ay 
basis, to preserve the symmetry of the walk. They obtain a spreading rate that is basi- 
cally classical, with a prefactor controlled by the rate of the measurements. The average 
squared distance from the start, measured by the variance, grows as a series of arcs at the 
quantum rate (quadratic), with regular resetting due to the measurements, pulling the 
overall rate back to linear (classical). They also generalise to random intervals between 
measurements, drawing analogy with Brownian motion. 

A unitary decoheren ce model affecting the position also appears in the same paper 



Romanelli et al.l (|200J) where they analyse a quantum walk in which links between the 



positions on the line are broken with probability p. The transition that should have taken 
place is turned into a self loop for that step of the walk. The unitary operator for each 
step is modified to take account of the links that happen to be broken at that time step. 
For low rates of link-breaking, the quantum behaviour persists, while for high rates the 
walk makes less progress than even an unimpeded classical random walk. 



5.4. Multiple coins in the walk on the line 



Another analytically tractable approach to reducing coherences in the walk on a line 
is to enlarge the size of the coin state space, and use parts of it in turn as the walk 
progresses. This has the advantage of the dyna mics remaining pure ly unitary, rendering 
the calculations simpler. The first such study, bv lBrun et al.l (j2003cf ). considered multiple 
coins used in sequence, with the sequence repeating after all the coins had been used once. 
This produces a quantum walk that still spreads linearly with the number of steps, but 
with the rate of spreading reduced inversely by the number of coins. Only if a new coin 
is used for every step of the walk does it become equivalent to the classical r andom walk . 



This i s in contrast with the behaviour obtained by decohering the coin, (IBrun et al 



2003a ). which always results in classical limiting behaviour, as pointed out in lBrun et al 



(|2003bfl . Classical behaviour is thus associated with an environment so large that one 
never comes close to the Poincare recurrence time over the timcscales considered. 

Rela ted studies that obs erve or exploit the behaviour when the coin space is limited 
include iFlitnev et al.l (j2004l ) who use multiple coins to create a Parrondo game, by having 
the amplitude of the coin flip depend on several previous coins rather than just one. 
Ermann et al.l (j2006l ) use an enlarged coin space to show explicitly that the behaviour 
changes from classical back to quantum once the size of the enviroi iment space is used 
up. They use the coins in a random order rather than sequentially as lBrun et al.l (|2003ar ) 
do. lRibeiro et al.l|2004 h study quasi-periodic sequences of coin operators using numerical 
simulation, and again find different spreading rates, but still an overall linear dependence 
on the number of steps. More realistic models of a finite-sized environment would have 
the environment degrees of freedom interacting with each other, but this has not yet 
been studied in the context of quantum walks. 
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Fig. 7. Numerical data for M^{p) on cycles of size = 48, = 49, iV = 50 and 
N = 51, for coin (dotted), position (dashed) and both (solid) subject to 
decoherence, using e = 0.002. Both axes logarithmic. 



5.5. Effects in the walk on the N-cycle 

Recall from W2.2\ that pure discrete-time quantum walks on cycles don't mix (except 
instantaneously for a few small special cases mentioned in i j2.3p . unless the time- averaged 
probability distribution, equation (jl3p is considered. In that case, mixing to the uniform 
distribution M(e), equation (fTS)) . does occur for some choices of coin operator, but this 
comes at the cost of requiring a number of steps linear in the inverse accuracy 1/e 
instead of l ogarithmic, unless a number of repetitions are combined in an amplification 



procedure. iKendon and Tregennal (120031 ) carried out numerical studies of decoherence 



in cycles. They evaluated M(e,p), the mixing time for P(x,T,p) c.f. equation (|15p . for 
walks on cycles of sizes up to ~ 80, both for pure states, and in the presence of the 
same types of decoherence as described in for the walk on a line. They did not 
include any amplification procedure. For odd-iV cycles with no decoherence, they report 
that M(e) ^ N/e as compared to the up per bound of M(e) ~ N\ogN/e^ given by 
Aharonov. D et al.l ( 200lh . iRichten ( 2006bl ) has recently confirmed this analytically. 

Under the action of a small amount of decoherence, the mixing time becomes shorter 
for all cases, typical results are shown in figure [7l If the coin operator is chosen such 
that the even-TV cycles do not mix to to the uniform distribution in the pure quantum 
walk, then the addition of decoherence causes them to mix to the uniform distribution. 
Although for N divisible by 4, the coin-decohered mixing time shows a minimum below 
the classical value at p ~ 2/7V, this mixing time is > N'^/32e, i. e., still quadratic in 
N. Thus, although noise on the coin causes the even-A^ cycle to mix to the uniform 
distribution, it does not produce a significant speed up over the classical random walk. 
Decoherence on the position produces a minimum mixing time Mmin(e,p) = 0{N/e), thus 
the even-A^ cycle mixes to uniform in linear time for a suitable choice of decoherence rate 
Pmin ~ 16/A^^, independent of e. 
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Fig. 8. Numerical data for M^{p) on cycles of size iV = 48 to iV = 51, for 
decoherence applied to coin only (dotted), position only (dashed) and both 
(solid), using e = 0.002. Both axes logarithmic. 



For all types of decoherence, the odd-A'' cycle shows a minimum mixing time at a 
decoherence rate somewhat earlier than the even-TV cycle, roughly p = but because 

of th e oscillatory nature of Pjx, T, p) , the exact behaviour is not a smooth function of p 
or e. ( Kendon and Tregennal ( 20031 ) illustrate this in their figure 4, we will show a similar 
example later.) As decoherence on the position (or both coin and position) increases, the 
oscillations in P(a;, T, p) are damped out. At p ~ 16 /N'^, the mixing time passes smoothly 
through an inflexion and from then on behaves in a quantitatively similar manner to the 
adjacent-sized even- TV cycles, including scaling as Minin(e,p) = 0{N/e) at the inflexion. 
Thus for < p < 16/iV^ there is a region where the mixing time stays linear in N. 

However, one can do better than this: decoherence causes the quantum walk to mix 
to the uniform distribution without any averaging of P{x,T,p), thus retaining the log- 
arithmic scaling with e for the mixing times. It is obvious that decoherence must do 
this, since high decoherence rates reproduce the classical random walk, which has this 
property. The only question is whether it does so effective ly enough to be useful, and 
the answer is t hat it does, for decoherence on the position ( Kendon and Malover . 200d : 



Richteii . l2006a ) . Figure [8] shows the mixing time M{e,p) corresponding to P{x,T,p) as 



a function of p for cycles of size 48 to 51. Since even-sized cycles only have support on 
half the positions at any one time, we have defined the uniform distribution only on 
the sites where it has support (we can easily fix this if necessary by averaging over two 
consecutive time steps). It thus behaves like an odd-sized cycle of half the size. Other 
than this technicality, the behaviour is the same, oscillating peaks end abruptly at the 
minimum mixing time, followed by a smooth rise to the classical value as the decoherence 
rate is turned up. The decoherence rate at the minimum mixing time is approximately 
7r/iV (odd-sized) or 2-k/N (even-sized), and the mixing time itself scales as 0(A^log(l/e)). 
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Fig. 9. Difference between the probability distribution P(x,T,p) and the uniform 
distribution expressed as the total variational distance, equation (fT3)) . for = 49 
for the case with decoherence on both coin and position. Both axes logarithmic. 
The value of e used in figures [7] and [8] is shown as a horizontal line. 



This thus provides a quadratic improvement over the classical mixing time, which scales 
as 0(Ar2iog(l/e)). 

If we take a look at the behaviour of | \P{x, t,p) — Pu\\tv, shown in figure [H for N = 49, 
we can see this too is fluctuating at the quantum end, with the period of fluctuation 
lengthening as it reaches the minimum. On the classical side of the minimum it increases 
smoothly to the classical value. The time at which the curves last cross the horizontal 
line at e = 0.002 is the mixing time as plotted in figure [H In the quantum regime, a 
different choice of e thus causes a jump in the value of Af(e,p) if it happens to touch 
the next peak in \\P{x,t,p) — Pu||tv This transition from underdamped (oscillating) to 
overdampe d (smooth) has been obs erved in a number of decoherence studies, such as the 
hypercube (jKendon and Tregennal 12003 1 , described in the next section, and the studies 
of decoherence in continuous-time quantum walks described in Sj6l The critical damping 
point also seems to be associ ated with the point as whic h all quantum correlations are 
destroyed by the decoherence ( Malover and Kendon . 2006h . just as this marks the optimal 
to p hat distributi on in the walk on the line. 



Richten (j2006b[ ) actually proves not only that the optimal mixing time is O ( log ( 1 /e) ) , 



but also that three diffe rent strategies for optimisi ng the mixing time are equally effective. 
The "warm start" from lAharonov. D et al.l ( 20011 1 can be viewed as running the quantum 
walk for several segments of random length separated by measurements. The decoherence 
applied with probability p is a quantum walk of length T interrupted by pT measurements 
at random times. Simply running the quantum walk for several segments of equal length 
with measurements applied at the end of each is also effective: the random outcomes 
from the measurements are sufficient to optimise the mixing. This shows that the effects 
of decoherence are quite robust, insensitive to the precise way in which it is applied. 
The overall conclusion is thus the same as for the walk on a line, there is a useful 
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window within which decoherence enhances rather than degrades the quantum features 
of the walk. Moreover, it is actually necessary to ensure the scaling with the precision 
e is efficient. Additionally, decoherence ensures the quantum walk mixes to a uniform 
distribution, regardless of the initial conditions or choice of coin operator. 



5.6. Effects in the walk on the hypercube 

Kempj (|2003bLl2005h picked up from where iMoore and Russel] |2002f ) left off and ana - 
ysed hitting times in discrete-time quantum walks on the hypercube. By hitting time we 
mean the time it takes for the quantum walk to reach the opposite corner from where 
it started from. Kempel considers two types of hitting times, one-shot, where a measure- 
ment is made after a predetermined number of steps, and concurrent, where the desired 
location is monitored continuously to see if the walker has arrived. While the one-shot 
hitting time is a pure quantum dynamics, the continuous monitoring of the concurrent 
hitting time measurement removes part of the quantum coh erences, thus it constitutes a 
type of decoherence. More recently, Krovi and BrunI ( 2006a ) have defined a hitting time 
in closer analogy with the usual classical hitting time, the average time of first arrival at 
the target location. 

Following Kempel . the one-shot hitting time -ffonc-shot(?'o) between \xq) and \xf) is 



shot(ro)-{r : |(a;/|U^|xo>P>ro} 



(78) 



where rg is a threshold probability of being at the target state, U is the evolution of one 
step of the walk, \xo) is the starting state and \xf) is the location being hit. 

To define t he re maining two hittin g times, we need to define a measured quantum walk, 
both Kempe . and Krovi and BrunI give essentially equivalent definitions. For hitting a 
single final location we need a measurement with two outcomes (this is easily generalised), 
say P/ and Qf — 1 — Pj, where P/ = \xf){xf \ (8) Ic, i-e., the projector onto the state 
for any coin state. If the walker is found at the walk is assumed to end, it acts 
as an absorbing boundary. If the walk does not reach \xf) after t steps, the state of the 
walk can be written 



(79) 



TriU^Qp] 

where Ugp = QfUpV^Q^ This differs from equation in that we keep only one 
outcome from the measurement and renormalise the density matrix. Since we know the 
measurement outcome, we know whether the walk has arrived or not, so we only need to 
keep the part that did not yet arrive at \xf). The first arrival (also called first crossing) 
probability after t steps can thus be written 



r{t) = Tr 

Kempel defines the concurrent hitting time 



He, 



iTcnt(?'o) = min{T \ 3 t < T : r{t) > tq} 



(80) 
(81) 
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The hitting time defined by iKrovi and BrunI is 



average 



t—oo 

E 

t=0 



trit), 



(82) 



the average arrival time. lKrovi and BrunI compare the concurrent hitting time with their 
average hitting time and find they both scale in a similar polynomial manner, compared 
with the exponentially long hitting time of a classical random walk, provided they use 
a Grover coin for the quantum walk, and consider \xq) — |0 . . . 0) and \xf) — |1 . . . 1). 
T heir nunierical results suggest the bound of ifconcurront(f^(l)) < O(n^logn^) obtained 
bv lKempel (|2003bll2005l) is not tight, the actual hitting time looks to be better than this 
in the large-n limit. 

Krovi and BrunI also investigate the effects of breaking the symmetry of the Grover 



coin, by using a DFT coin, and by distorting the hypercube lattice. The DFT coin has the 
property that for certain configurations the hitting time becomes infinite, while a Grover 
coin with a small distortion of the hypercube lattice increases the hitting time somewhat, 
though it is still less than the classical random walk hitting time. A nalogous results 
for co ined quantum wa lks on the "glue d trees" graph were noted by iTregenna et al 



(l2003l) . Further work bv lXrovi and BrunI |2006b |) on the question of whether and under 



what conditions quantum walks show dramatically different properties (speed up or slow 
down) compared with classical random walks suggests that it is highly dependent on the 
symmetry of the g raph. For the Grover coin, this is exemplified by the quantum walk 
search algorithm of Shenvi et al. ( 20031 ) (described in ^4.1^ . where any disturbance of the 
symmetry causes the walker to converge on the marked state. 



Kendon and Tregennal ([20031) investigated the effects of adding extra decoherence to 



the quantum walk on a hypercube over and above that implied by the measurements of 
the target location. For both one-shot and concurrent hitting times (this work predates 
Krovi and Brun ( 2006ah ). the key parameter is the probability r{t) of finding the walker 
at the chosen location. Their numerical calculations show that all forms of decoherence 
have a similar effect on r{t), see figure [TUl reducing the peaks and smoothing out the 
troughs. For the one-shot hitting time this is useful, raising r{t) in the trough to well 
above the classical value, so it is no longer necessary to know exactly when to measure. 
For p < 1/n, the height of the first peak scales as rp{t) — r{t){0) exp{— (n -|- a)p}, where 
^ Q; ^ 2 depending on whether coin, position or both are subject to decoherence. 
An exponential decrease in the presence of decoherence sounds about as bad as it could 
reasonably be, and for long times, of course, decoherence reduces the walk to classical 
behaviour. However, the hitting times are short, only ~ m:/2 steps, so p ~ 1/n only 
lowers r(t) by a factor of 1/e. For algorithmic purposes this is insignificant, only a factor 
of order unity and thus still exponentially better than classical. Note also that the size 
of the graph (measured in number of nodes) is 2", so the decoherence has only a linear 
effect measured in terms of the size of the graph. 

The concurrent hitting time already includes a portion of decoherence, no extra features 
are produced by the addition of unselective decoherence, but there is still a range of 
< p ^ 1/n within which the quantum speed up is preserved. Note that in both the 
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time step 



Fig. 10. Hitting probability on a 9-dimensional hypercube for one-shot (left) and 
concurrent (right), perfect walk (circles), with p = 0.05 (dotted), p = 0.1 ~ 1/9 
(solid). Classical hitting probability barely visible (dashed). 



one-shot and concurrent cases, p ~ 1/n is a critical damping rate, smoothing out the 
second peak (shown at around 40 (= 3n7r/2) steps in figure fTO)). 

Decoherence in discrete-time walks has thus provided us with a number of notable 
common features, such as smoother spreading, enhanced mixing, and a transition from 
under-damped (quantum) to over-damped (classical) behaviour analogous to classical 
damping in a harmonic oscillator. We will see in the next section that these features also 
appear for the continuous-time quantum walk under decoherence. 



6. Decoherence in continuous-time walks 

Less work has been done on decoherence in continuous-time quantum walks, probably 
in part because the numerical simulations require more resources (to integrate rather 
than iterate the dynamics). Analytical calculations for the continuous-time walk are 
often simpler than for the discrete-time walk due to the lack of a coin, but in the case 
of decoherence this removes the simple case of decoherence on the coin only that has 
been so successfully studi ed in the discret e -time walk. Nonetheless, there are two notable 



recent analytical results. iFedichkin et al.l (|2006l ) obtain the scaling for mixing times for 



the walk on the cycle for both small and large decoherence rates. They augment this 
with numerical studies showing a minimum mixing time for intermediate decoherence 
Alagic and Russell ( 2005h solve the decoherent walk on the hypercube for all rates 



rates. 

of decoherence. These results will be discussed in separate sections below. First, we will 
set up a general model for decoherence in continuous-time quantum walks analogous to 
the basic model used for discrete-time walks. Then wc will briefly examine the effects of 
decoherence in the continuous- time quantum walk on the line, obtained numerically, to 
illustrate the key features we can expect to see. 

The pure state evolution of the continuous-time quantum walk given by equation (|28p 
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can be expressed in density matrix form as 

dpjt) ^ 
dt 



T-l [A, p] , 



(83) 



where [a, h] = ah — ha denotes the commutator. We can add a non-unitary decoherence 
to this in exact analogy to the discrete-time method by writing 



p{t + St) = (1 - pSt) {p{t) - tjSt [A, p]) + pStVp, 



(84) 



where V represents the action of the noise, for example, a set of projectors {Pj}, such 



that Vp = J2, 



^jpPj- For uncorrelated noise, taking St 



this yields 



dpjt) 
dt 



(85) 



The effect of the extra two terms is to reduce some or all of the off-diagonal elements of 
p{t) at a rate p per unit time. 



6.1. Effects in the walk on the line 

In contrast to the extensive studies of the discrete-time walk on the line, there are no 
published analytical treatments of decoherence in a continuous-time quantum walk on 
the line. We can easily investigate numerically: for uncorrelated noise events at a rate p 
per unit time, we can unpack equation (|85p into 

d 



dt 



Px,yit) 



Px,y + 1 — Px + l,y — Px-l,y + Px,y-1 



P (1 - Sx,y) Px,', 



(86) 



where x,y £ Z . ( Fedichkin et al. ( 2006 ) credit this model as developed by Gurvit j (1997); 
Gurvitz et al.l (j2003l ).) Figure [Tl] shows the results of evolving a continuous-time quantum 
walk on the line with hopping rate 7 = 0.5 until t = 40, for various decoherence rates. 
This is to be compared with figure [l] which shows the same type of study for a discrete- 
time quantum walk. The continuous-time walk does not produce quite such a nice top hat 
as the discrete-time walk, partly due to having support on all sites rather than alternate 
sites, though it would probably be equally useful in practice. There is, however, a striking 
new effect: unlike with the discrete-time walk, where decoherence rates larger than one 
have no useful interpretation, there is no reason why p in equations (j85|) and (|86p cannot 
be made arbitrarily large. In figure [HJ p ~ 8000 is shown: the walk hardly manages to 
leave the starting point. High values of p can be interpreted as continuous monitoring, 
which projects the walk into the initi al state with high probabilit y. In other words, the 
walk suffers a quantum Zeno effect ( Misra and SudarshanL 1977 ). a phenomenon that 
does not appear in coined quantum walks. One consequence of this is that quantities 
such as the concurrent hitting time are infinite, since the continuous monitoring will 
ensure the continuous-time quantum walk never arrives. 



6.2. Effects in the walk on the N-cycle 

The time-dependent non-unitary evolution of p( t) on a cycle is also gi ven by equation 
(|86p. only with x,y £ for a cycle of size N. Fedichkin et al. ( 2006h showed how to 
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-20 20 

position on line 



Fig. 11. Probability distribution for a continuous-time quantum walk on the line 
with 7 = 0.5, at time t = 40, for various decoherence rates p give in the key. 



solve this for small and large values of decoherence. We will now outline their method. 
The following substitution 

R..y = (87) 
converts the differential equations into a set with real coefficients, 

-^Rx,y — ^ (Rx.y+l + Rx+l,y — Rx-l,y — Rx,y-l) " P (1 " ^x,y) Rx,y (88) 

Note that this leav es the diagonal e lemen ts of p identical to the diagonal elements of 
R. As a side note, iFedichkin et al.l (j2006f ) observe from these equations that if p = 0, 
there is an exact mapping of the quantum walk on a cycle onto a classical random 
walk on a two-dimensional torus, and if p ^ 0, there is still an exact mapping of the 
quantum walk on a cycle onto some classical dynamics on a directed toric graph. Similar 
mappings (one-dimensional quantum to two-dime nsional classica l) also occur in quantum 
phase transitions in spin chains (see, for example, Sachdev . 1999[ l. and can prove a useful 
technique for analysing quantum systems, though we will not need to use it explicitly 
here. 

Equation (j88p can be expressed as a linear operator equation 

^R(t) ^ nL+pV) RU). (89) 
at 

One can think of L and U as iV^ x iV^ matrices and R{t) as a vector of size N'^ (instead 
of a matrix of size N x N). The indices of R{t) as a vector become xN + y where x, y 
are the pair of indices of R{t) as a matrix, but we will write xy for such indices to 
reduce complexity in the notation. The formal solution of ([89l) can be written R{t) = 
gt(L+pU)^|'Q-j^ The explicit solution for arbitrary decoherence rates has not been obtained, 
but one can examine the behaviour of the mixing time under the action of small and large 
decoherence, using standard perturbation theory methods. When the decoherence rate p 
is small such that pN <C 1, pV is treated as a perturbation of L. The unperturbed linear 
operator L (which evolves the pure state quantum walk in density matrix formalism) has 



V. Kendon 



38 



eigenvalues 



\nn — * sin 



7r(m + n) 
N 



cos 



77(171 — n) 
N 



with corresponding eigenvectors 



Vf,u,mn = -f-f exp —{mfj, + nv) 



N 



2TTi 



N 



(90) 



(91) 



Fedichkin et al. 2006. for 



Here /lu, mn are indices < 11, v, 771,71 < N combined like x, y above. With due con sid- 
eration of the degeneracies in the eigenvalues Xmn of L (see 
details), the eigenvalue perturbation A„j„ turns out to be 

X J^-2) 



N 



Thus, the solution is of the form 



exp 



27ri 

— [mx + ny) 



(92) 



(93) 



The probability distribution of the continuous-time quantum walk is given by the diagonal 
terms (considering R{t) a,s a.n N x N matrix again) P{x,t,p) = Rx,x{t), that is 



P{x,t,p) = ^ + ^ X! ^ V+"](mod Ar),o) X 



(m,n) 



X exp 



it sin 



7r(TO + n) 
N 



cos 



m,n 

7r(m — n) 
N 



exp 



+ (1 - '5m,«)e-P" 
27ri 



N 



(m + n)x 



(94) 



Fedichkin et al. ( 2006h then calculate an upper bound on the mixing time M{e). Defining 

(95) 



N-l 



where lon — exp(27ri/A^), and noting that 



N-l 



N-l 



■mx 
N 



(96) 



m.n— m— 

and that |Al2,(t)| < 1, we can simplify equation (|94p as follows. 



P{x,t,p) 



1 

TV 



Ml{t/2) 
1 



-tp/N _ I 



N 



N 

(1 - 2/7V) 



M2x{t) 



2-{N mod 2) 



N 



(97) 



Adding a now-trivial summation over x gives us the total variational distance required 
to obtain the mixing time: 



N-l 

E 

x=0 



< e-P^' (iV + e-*P/^-l). 



(98) 
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Since e < 1, the above equation shows that the mixing condition is A^e ^ « * < £■ 
This gives the mixing time bound of 

Mie) < ^log(f)[l + ^]. (99) 

Some comments on this result are in order. First of all, this is the mixing time for the 
instantaneous probability distribution, so it should be compared with the discrete-time 
results in figure [51 Just as for the discrete-time walk, the first key effect of decoherence 
is to cause the continuous-time quantum walk to mix. On the face of it, M(e) appears to 
scale as log(A^), which would imply an exponential speed up over classical mixing times 
of 0{N'^). However, this formula is only valid for pN ^ 1, the 1/p scaling thus ensures 
that M{t) ^ N for the range of validity of the result. What this result tells us is that 
decoherence causes the continuous-time quantum walk to mix with the dominant effect 
a scaling of and only weak (logarithmic) dependence on N and e. It does not tell us 
anything about the optimal decoherence rate that might give a minimum mixing time, 
nor what value that minimum mixing time might take. 

We can calculate the time-averaged mixing time from the results for P{x,t,p) above, 
and compare it with the discrete-time and classical values. Again, we will obtain a result 
only valid in restricted ranges of parameters. Working with equation (|97p . and assuming 
iV 3> 1 for simplicity, we obtain 



P{x,T,p)- 



< 



1 



-pt 



1 

N 



dt 



1 

N 



< 



pT 



(1 



(100) 



Summing over x to obtain the total variational distance trivially multiplies by N. Com- 
paring this to e as before gives 



-d 

pT^ 



-pT 



)<e 



(101) 



Assuming pT ^ 1 (which is consistent with our other assumptions, 3> 1 and pN <ti 1, 
it requires T 2> iV) gives 



N 

M{e) > -, 
pe 



(102) 



which is similar to the scaling found numerically for the discrete-time quantum walk. 
In particular, the scaling with N and e is no longer logarithmic, as a consequence of 
the time-averaging. The scaling with p is still inverse, indicating that, in general, pure 
continuous-time quantum walks on cycles do not mix to within e of uniform, even when 
time-averaged distributions are considered, as already noted in ijS.SI 

We now turn to the analogous calculation for large decoherence rates. When the deco- 
herence rate p is large, that is, when p ^ 1, the quantum walk is experiencing a quantum 
Zeno effect that suppresses most of the transitions. We can therefore neglect all but the 
first off-diagonal terms in the matrix R{t) from equation (|88p . Retaining only matrix 
elements that are of order 1/p produces a truncated set of differential equations for the 
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elements along the major and the two adjacent minor diagonals: 



1 



'x,x+l 



{Rx+l,x+l ~ Rx,x) ~ pRx,x+l: 
{Rx,x -^x— — l) pRx^x — l' 



^x,x — l ^ 

To facilitate our subsequent analysis, we define 

Ojx — Rx^xt dx — Rx,x+1 ~t~ Rx+l,X' 

Then, we observe that 

/ _ {dx — dx-i) ,/ _ {cix+i — o-x) _ , 

The solution of differential equation at large p has the form 



Af-l 



sm 



exp 



2 Try 
N 



j/=0 



2p 



(103) 
(104) 

(105) 

(106) 
(107) 

(108) 



Based on the above analysis, the full solution for R{t) is given by 

Ox if X = 2/ 
Rx,y{t) = { dx/2 if \x -y\ = l 
otherwise 



(109) 



The total variation distance between the uniform distribution and the probability 
distribution of the decoherent quantum walk is given by 



ax{t) - 



JV-l 

E 

x=0 

which simplifies to 

E 



JV-l 

E 

x=0 



Af-l 



ivE«-p 



sm' 



2 Try 
N 



V=0 



2p 



t exp 



/ 2Trixy 



1 

N 



x=0 



axit) - 



N~l 

N E 



x=0 



N-1 

E 

y=l 



2 Try 



2p 



t cos 



/ 2Tryx 



(110) 



(HI) 



A lower bound on the mixing time for large decoherence rate p can be derived as 
follows. Note that 



N-l 

E 

x=0 



ax{t) 



1 

N 



> 



ao{t) 



1 

iV 



JV-l 



N 



^exp 



2 Try 
^ N 



2p 



> ^exp 



2p 



t 



(112) 



(113) 



where the first inequality uses the term .t = only and the second inequality uses the 
terms y = 1, iV — 1. This expression is monotone in t, and is a lower bound on the total 
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variation distance. It reaches e at time Tiower, when 

Slower 



Sin' 



2p / 2 



2pN^ ( 2 



(114) 



for large ^ 1. 

An upp er bound on the mixin g: time for large decoherence rate -p can also be derived 
(details in lPedichkin et all |2006h ) 



■In 



2 + e 



(115) 



(116) 



For large decoherence rates p 3> 1, the mixing times are bounded as 

These bounds show that M (e) is linearly proportional to the decoherence rate p, but is 
quadratically dependent on iV, and logarithmically dependent on e. 

The change in the form of the scaling with p from reciprocal for small p to linear for 
large y suggests there is a minimum mixing time at some intermediate decoherence rate. 
Fedichkin et al.l (j2006l ) confirmed this numerically, reproduced in figure fT2l 

The effects of decoherence on continuous-time quantum walks on cycles is thus similar 
to the e f fects on discrete-time walks on cycles, as described in ^ 35. 51 In fact, the proofs in 
Richteii (j2006b[ ) are sufficiently general they apply equally to the continuous- time walk as 
to the discrete-time walk. This fills in the gap between small and large decoherence rates, 
and proves that continuous-time walks, with 0(log(l/e)) decoherence or measurement 
events, mix in time 0(iVlog(l/e)) on the cycle (and the d-dimensional torus). 
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6.3. Effects in the walk on the hypercube 

The hypercube is an attractive subject for analytical calculations of quantum walks 
because, under certain conditions, it offers two routes to simplify the analysis. If the 
walk has the appropriate symmetry, it is constrained to a subspace that maps to a walk 
on a line segment with positions corresponding to the Hamming weight of the node labels 
(with position dependent transition rates). Alternatively, the continuous-time walk can 
be factored in Fourier space into independent walks on each qubit of the labels of the 
locations, as explained in ^3.41 Furthermore, both mixing properties for spreading over the 
whole graph, and hitting properties for traversing to opposite corners are non-trivial on 
the hypercube. We thus have a large number of comparisons to make between discrete and 
continuous-time quantum walks, and classical random walks. lAlagic and Russelll ^200^ 
provide a complete solution to the dynamics of the continuous-time quantum walk on the 
hypercube subject to decoherence. We will sketch their method of solution, then discuss 
t heir results. 

Alagic and Russelll give their quantum walk an energy k - this is equivalent an arbitrary 
hopping rate j = k/n instead of the choice of 7 = 1/n for a graph with vertices of 
degree n. We consider the continuous quantum walk on the n-dimensional hypercube 
with energy k and decoherence rate p, starting from the initial wave function = 
|0)®", corresponding to the corner with Hamming weight zero. The decoherence operators 
in equation ([85|) project the walker onto one of the vertices of the hypercube chosen 
uniformly at random. It is useful to write these projectors explicitly in terms of single 
qubit projectors Hq and Hi onto |0) and respectively. We have 

„ 1 ' 



l<j<n 



Hi 



>m 



(117) 



where Hq = 1 (g) ■ ■ ■ 1^ 1 (E)IIq (g) 1 (g) ■ ■ ■ (g) 1 with the non-identity projector appearing in 
the jth place and similarly for H{ . 

We now show that, with the model of decoherence described above, each dimension 
still behaves independently. Recall equation ([1T|) showing that the adjacency matrix for 
the hypercube decomposes into a sum of tensor products each acting only on a single 
qubit. Since the noise operators Pj also have this structure, provided p(0) starts off in a 
state that is also decomposable (such as the all-zero vertex) , the subsequent evolution will 
maintain this structure and will remain decomposable as a system of n non-interacting 
qubits. Since the qubits are now in mixed states Pj, the superoperator acting on them is 
equivalent to a 4 x 4 matrix 



J 
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p(l g)t)+ p(Hi «) Hi) + p(Ho g) Ho)] 



(118) 



gjto-a^j which is exactly the superoperator 
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formulation of the dynamics of the non-decohering walk. So far we have shown that the 
walk with decoherence is still equivalent to n non-interacting single-qubit systems. We 
now analyse the behaviour of a single-qubit system under the superoperator e-^. The 
structure of this single particle walk will allow us to then immediately draw conclusions 
about the entire system. 

-pt-ipt -pt+ipt 



The eigenvalues of J are 0, — ^ 



Here /3 = y/l6k'^ — is a 



2n 2n 

complex constant that will later turn out to be important in determining the behaviour 
of the system as a function of the rate of decoherence p and the energy k. The ma- 
trix exponential of in this spectral basis can be computed by inspection. To see how 
our superoperator acts on a density matrix po, we may change po to the spectral ba- 
sis, apply the diagonal superoperator to yield p(t), and finally change p{t) back to the 
computational basis. At that point we can apply the usual projectors Ho and Hi to de- 
termine the probabilities of measuring or 1 as a function of time. The result (details in 
Alagic and Russell . 2005h is 



r n 1 1 f 

P[0] ^ - + -e-- 



0t\ p . /(3t 
cos \ — ] + — sm — 
2nJ (3 \2n 



P[l] 



cos 



+ 



/3 



(119) 



All the properties of the decohering walk can be deduced from these equations. Just as 
was observed in the mixing behaviour of cycles, the decoherence exhibits three distinctive 
regimes analogous to underdamping, critical damping and overdamping in a damped 
oscillator. We will now describe the properties of these regimes. 

The underdamping regime occurs for p < 4fc. To determine the mixing times, we solve 
P[0] = P[l] = i giving exact instantaneous mixing times at 



n(27rc — arccos(p^/8fc^ — 1)) 



yieP 



(120) 



for all c G Z, c > 0. At these times, the total variational distance between the walk 
distribution and the uniform distribution is ze ro, but the walk "u n mixes " again as time 
progresses. This corresponds to the result in iMoore and Russelll (|2002r ) for pure state 
quantum walks, and extends it to regions with small decoherence rates. These mixing 
times scaling linearly with n represent a quantum improvement over the classical random 
walk mixing time of Q(nlogn), with the caveat that the classical mixing time is not 
instantaneous, so strictly speaking we are not comparing like with like, in the quantum 
case one must sample at exactly the right moment in time. Note that these periodic 
mixing times decay with p and disappear altogether when p > 4fc, so the qu antum 
advantage tails off before the critical damping point is reached. Richter ( 2006bh offers 
a weak bound of 0(71"^/^ log(l/e)) for the mixing time at around the critical damping 
point, which is slower than classical. 

By examining the local maxima of P[l], we can determine that the walk has approxi- 
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mate instantaneous hitting times to the opposite comer (1, 

27rn(2c+ 1) 



, 1) at times 



one-shot 



(r) 



(121) 



■y/l6fc^ — 

for aU c G Z, c > 0. However, the probabihty r(p, c) of measuring an exact hit decays 
exponentially, 

'1 1 



r{p, c) = 



:-shot 



(a 



K7r(2c+1) 



(122) 



and it does 



When no decoherence is present, the walk hits at iJonc-snoi — ni. 

so exactly, i.e. r(c) = 1. This corresponds to the result of Kempel ( 2003b . 2005 ) for 
the discrete-time quantum walk, see H5.61 thus both discrete-time and continuous-time 



quantum walks on the hypercube exhibit exponentially faster hitting times than classical 
random walks. The difference is that the continuous-time quantum walk hits exactly, 
while the discrete-tim e quantum walk hits with a probability that is less than one, see 
Kempel |2003bl boOSh for details. 



For a fixed p > 4fc, lAlagic and Russell (120051 ) show that the walk behaves much like 
the classical walk on the hypercube, the measurement distribution of the walk converges 
to the uniform distribution in time A/(e) ~ 0(nlogn), just as in the classical case. 

As p — > cx), the walk suffers from the quantum Zeno effect. Informally stated, the 
rate of decoherence is so large that the walk is continuously being reset to the initial 
wave function |0)®" by measurement. By inspection of equations (|119p . it is clear that 
P[0] -> 1 and P[l] ^ as p ^ cx). 

This concludes our overview of decoherence in continuous-time walks. We have ob- 
served some of the same effects of enhanced mixing and smoother spreading, and also 
shown that where the decoherence doesn't directly help, the quantum speed up is at least 
robust under small amounts of decoherence. 



7. Quantum vi^alks in physical systems 



As well as the potential for algorithmic applications that has been in the background of 
most of our discussion of decoherence in quantum walks in this review, physical systems 
can be made to perform a quantum walk, and this process is of interest in its own right 
as an example of precise coherent control. Indeed, the first suggested application of a 
quantum walk in a physical syst em was to test decoherence rates in a single trapped ion. 
Travaglione and MilburnI |2002h describe how a coined quantum walk can be performed 
using the vibration modes of the trapped ion coupled to an internal state of the ion for 
the quantum coin. They note that a cycle of size j V = 4 performed in this way should 



exhibit periodic behaviour, see §2.31 1 Sanders et al.l (j2003l ) suggest a quantum walk using 
the phase of the electromagnetic field in a cavity as the position (so the walk is again on 
a cycle) with the role of the coin played by an atom, also in the cavity, which couples 
to the cavity field depending on its internal state. Since the cavity field cannot be in a 
state localised at a point i n phase space, the positions must be approximated by small 
regions that do not overlap. ISanders et al.l provide a detailed analysis of the ex perimental 
parameters necessary to maintain distinct positions during the quantum walk. iDiir et al 
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( 20021) jropose a quantum walk on the line using an atom hopping between traps in an 
optical lattice. The role of the coin is again played by the internal state of the atom. 
One of the potential sources of error in this system is imperfect transitions between 
lattice site s, so, along with analysing errors due to imperfect Hadamard operations, 



Dtir et al.l also consider the effects of imperfect shift operations. Their work is numerical. 



but shows that imperfect shifts, or spontaneous tunnelling, smear out the distribution 
without necessarily reducing the quantum features. This analysis is typical of the extra 
factors that need to be considered in experimental systems. Since the quantum walk 
takes place in a subspace of the full quantum system, many of the potential sources 
of decoherence take the system outside of the subspace of the quantum walk, so the 
resulting decohered dynamics is not a classical random walk. All of these proposals are 
within current experimental capabilities but, to date, none has bee n carried out. 

The above examples are all of discrete-time coined quantum walks. lSolenov and Fedichkin 
|2006al lbl) analyse a ring of quantum dots in which the electrons in the dots perform a 
continuous-time quantum walk. They apply similar perturbation theory techniques to 
analyse the effects of decoherence as were described in §6.21 

There are a number of ways in which photons can be made to follow a quantum 
walk dynamics, and while some have been implemented, all have used high light in- 
tensities that admit a classical wave description. The earliest example, wh i ch pr edates 
the algorithmic discussions of quantum walks, is from iBouwmeester et al.l (|l999l ) , who 
present th e twin peaked distri bution of a quantum walk on the line without identifying 
it as such. Knight et al. ( 2003h discuss t his experiment, arguing that quantum walks on 
the line are actually not q uantum at all. iKendon and Sandera (|2004l ) explain how what 
Bouwmeester et al.l (|l999l) have done can be regarded as a quantum walk and suggest 



an enhanced experiment, in which monitoring the individual photon "walkers" reduces 
the dynamics to a classical random walk, to distinguish a classical "wave walk" from a 
quantum walk. 

A quantum walk experiment carried out by iRvan et al.l (j2005r ) uses a 3 qubit NMR 
system to perform a quantum walk on a cycle of size TV = 4. This is actually a quantum 
computation of a quantum walk, since the three qubits are used to represent the binary 
number labelling the vertex (two qubits), and t he qubit coin. T his is distinctly different 
from the p hysical syste ms described above, see Kendon ( 2006bh for more discussion on 
this point. iRvan et al.l also apply decoherence artificially to the qubits, to observe that 
the quantum walk degrades into the classical walk. 



8. Quantum walks on directed graphs 

Providing a quantum walk dynamics on undirected graphs is straightforward since the 
adjacency matrix of an undirected graph is symmetric, hence it can be turned into a uni- 
tary (with added coin) or Hermitian operator for any such graph. In contrast, determining 
whether a directed graph allows a unitary dynamics that re spects the graph structure 
is a hard problem, to which a general answer is not known. lAharonov. D et al.l (|200lh 
give a useful definition of how a quantum walk dynamics should re spect the st r uctur e 
of a graph, though other approaches are possible, see, for example, Montanajol ( 2007 ). 
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Severini ( 20031 2006 ) provides some necessary conditions, by studying the pattern of a 



unitary matrix, defined as 



1 if U,. 



^ and = if U^j 



= O.UA = U^P^ for 



some [/, then a quantum walk can naturally be applied to the graph with adjacency ma- 
trix A. The situation is simpler for some restricted types of graphs, such as line graphs. A 
line graph L'-"-^-' is obtained from graph G by making each edge in G a vertex in L(<=), and 
edges for each pair of edges in G that are adjacent. If x, y, z are vertices in G connected 
by directed edges {xy) and {yz) which can be traversed only from 2: to y to z, then the 
vertices in L^'^^ correspond to all {xy) S G, and the edges in L^*^^ correspond to all 
{{xy)(iiz)) where {xy) and {yz) are edges in G. This definition can be ap plied recursively 



to gen erate a family of line graphs from a single initial graph G, and iPakohski et al 
(|2003 h provide conditions under which line graphs can be quantised. 



Starting with an undirected graph , a corresponding directed graph can be constructed 
by a mapping due to lWatroua (|200l[ ). similar to the line graph definition above, only each 
undirected edge is treated as a pair of directed edges. The resulting directed graph is 
always the pattern of a unitary matrix. The corresponding dynamics, using an associated 
unitary matrix as the step operator, corresponds exactly to the coined quantu m walk 
dynamics (with a Grover coin) on the original undirected graph |Watrousl . [200i ). 

Directed graphs t hat are not the pat tern of a unitary matrix may still admit a quan- 
tum walk dynamics. iMontanard (j2007l ) defines "reversibility" for a directed graph, and 
proves this allows for a unitary quantum walk in which the coin selects between different 
sets of cycles in the graph, rather than different paths from each vertex. For irreversible 
graphs, Montanaro i suggests splitting the graph into subgraphs that are reversible, on 
which the quantum walk dynamics can be unitary, combined with a non-unitary step (a 
partial measurement) to switch between subgraphs in an irreversible way that nonethe- 
less preserves coherence within each reversible subgraph. If one is willing to all ow some 
transiti ons that are between vertices not connected by edges, the method of ISzegedv 
(|2004a[) also provides a quantum Markov process on directed graphs. 

This brief section only touches on the range of problems in the realm of directed graphs, 
reflecting the limited work done in this area. For example, no study of decoherence in 
quantum walks on directed graphs has been done, nor consideration of the class of graphs 
that might be the pattern of a CP map (rather than of a unitary matrix). 



9. Summary and outlook 

We have seen in our broad-ranging review of decoherence in quantum walks that the ef- 
fects are both non-trivial and have potential algorithmic applications. As the decoherence 
rate is increased, for moderate timescales the quantum features persist, usually to a point 
of "critical damping" after which the behaviour resembles that of a classical system. In 
the quantum region before the critical damping point, the useful algorithmic properties, 
such as spreading and mixing, are often better than the pure quantum dynamics. Indeed, 
for cycles and tori, some element of non-unitary evolution is essential for efficient mixing 
behaviour. Our tentative conclusion is thus that decoherence is a useful enhancement, 
and non-unitary dynamics are a promising avenue of enquiry in the quest for practi- 
cal applications for quantum computing. The incomplete nature of our conclusions are 
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in part a reflection of our limited knowledge of the conditions under which a quantum 
speed up can be obtained. 

The idea that the addition of classical randomness can enhance a useful algorithmic 
property also occurs in other settings, both qu antum and classical . An example from th e 
field of quantum information may be found in lKraus et al.l l|2005l ): lRenner et all |2005h . 
where the security of a quantum secret key is increased by adding a small amount of 
randomness to the key during the privacy amplification stage of the protocol. This can be 
understood in two ways. The goal is to reduce the amount of information available to an 
eavesdropper to a trivial level. The added randomness reduces the final size of the secret 
key, but it reduces even more the information that the eavesdropper has. This works for 
both classical and quantum keys. From a quantum perspective, the randomness takes the 
form of entanglement with an environment that the eavesdropper does not have access 
to. Since the total amount of entanglement is limited, this reduces the entanglement 
of the key with the eavesdropper, and thus reduces the information available to the 
eavesdropper. The interpretation of this in the context of quantum walks is of course 
rather different, but the common theme is that the decoherence or randomness induces a 
more uniform distribution, in this case in the form of a smoother "top-hat" distribution 
in the walk on the line, and faster mixing times on finite graphs. Some further insight 
is provided in iKendon and Maloveii (|2006r ). where the entanglement between the coin 
and the position is used to gauge the impact of the added randomness. The optimal 
decoherence rate turns out to just remove all the quantum entanglement, and can be 
thought of as removing the "extra" quantum correlations that would otherwise keep the 
distribution away from uniform. 

As will have become clear to the reader by now, there are many gaps in the analysis of 
decoherence in quantum walks, and, more importantly, many gaps in our understanding 
of what might make a useful quantum walk algorithm. Some of these shortcomings will 
only be accessible to serious study when (if) we have a larger (that is, larger than we 
can simulate classically) quantum computer available to test the performance of actual 
algorithms. In the meantime, there is still much useful progress that can be made, in 
particular in understanding the role of symmetry in the properties of quantum walks. 
Many of the examples studied so far turn out on closer inspection to be exceptional cases, 
a slight change in the symmetry of the dynamics or initial state can reverse the interesting 
quantum properties. The role of decoherence as a "memory eraser" to eliminate the 
dependence on the initial state, while still retaining the interesting quantum properties, 
is thus one of its most appealing features. My hope in offering this review is that it will 
be a useful summary of the current state of research in a field that has seen significant 
recent activity and advances, and that it will facilitate and spur further efforts in this 
area. 

I gratefully thank my colleagues for many useful and stimulating discussions of quan- 
tum walks, especially Dorit Aharonov, Andris Ambainis, Todd Brun, Ivens Carncirq^, Hi- 
lary Carteret, Andrew Childs, Richard Cleve, Jochen Endrejat, Ed Farhi, Will Flanagan, 
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